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PREFACE 


This  is  one  of  three  volumes  of  the  final  report  on  "Aircraft  Configuration  Noise 
Reduction”  submitted  by  the  Boeing  Commercial  Airplane  Company,  Seattle, 
Washington,  98124,  in  fulfillment  of  Department  of  Transportation  contract 
DOT-FA74WA-3497,  dated  1 August  1974.  This  work  was  completed  for  the  ATC 
Airport  Facilities  section  of  the  Federal  Aviation  Administration  (DOT).  Mr.  H.  C.  True 
was  the  Contract  Technical  Monitor. 

The  report  is  divided  into  three  volumes  for  easy  use  as  shown  below: 


Vol.  Ill 
computer 
program 
source 
listing 


• FORTRAN 
source  code 


Vol.  II 
computer 
program 
user's  guide 
and  other 
appendices 


• A-Computer  program  user’s  guide 

• B-Flight  effects  test  plan 

• C.  ..H 


Vol.  I 

engineering 

analysis 


• Theory 

• Test  review  and  data  analysis 

• Estimation  procedure 

• Conclusions  and  recommendations 


This  report  is  volume  H of  the  series  and  was  prepared  jointly  by  the  Noise  Technology 
Staff  of  the  Boeing  Commercial  Airplane  Company  and  the  Noise  Systems  Group  of 
Boeing  Computer  Services,  Inc.  This  volume  contains^ 

• **The  user’s  guide  for  the  computer  software  of.  the  '^Aircraft  Configuration  Noise 

Reduction"*8tudy. 

• ^A  preliminary  test  plan  for  assessing  forward  velocity  effects  on  wing  and  fuselage 

shielding.  , 

• Various  curves,  derivations,  and  background  theory  in  support  of  material 
presented  in  volume). 
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APPENDIX  A 


FAA-AIRCRAFT  NOISE  SOURCE  COMPUTER  PROGRAM 

USER’S  GUIDE 

D.  G.  Dunn,  D.  J.  Cecil,  and  L.  M.  Butzel 

1.0  SUMMARY 


This  appendix  is  intended  to  serve  as  a user’s  guide  to  the  integrated  system  of 
computer  programs  developed  in  accordance  with  the  requirements  of  contract 
DOT-FA74WA-3497.  The  computer  software  for  the  estimation  procedures  of  flyover 
noise  time  histories  developed  under  contract  NAS2-6969  has  been  used  where 
compatible,  thus,  retaining  the  consistency  in  the  input  description  found  in 
NASA  CR1 14650,  "Aircraft  Noise  Source  and  Contour  Computer  Programs  User’s 
Guide.”  Moreover,  since  the  computer  program  is  an  extension  of  the  program  developed 
under  contract  NAS2-6969,  this  appendix  may  be  considered  as  an  updated  version  of 
sections  1.0,  2.0,  and  3.0  of  NASA  CR1 14650,  structured  within  the  requirements  of 
contract  DOT-FA74WA-3497.  For  convenient  use  the  additions  necessary  to  incorporate 
the  effects  associated  with  the  shielding  of  noise  sources  near  a wing  surface  are 
indicated. 


Illustrations  provided  in  the  form  of  macro  flow  charts  show  the  computer  software 
developed  under  DOT-FA74WA-3497  and  its  linkage  to  that  of  contract  NAS2-6969. 
Sample  cases  of  input  data  are  included  with  a limited  sample  output  to  assist  the  user 
in  the  use  of  the  program.  The  programs  are  stored  on  magnetic  tape  and  are  written  in 
FORTRAN  IV  for  the  FTN  compiler  to  run  on  a CDC  or  CYBER  computer,  operating 
under  the  KRONOS  2.1  operating  system. 
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2.0  INTRODUCTION 


The  computer  program  developed  under  contract  DOT-FA74WA-3497  simulates  a 
procedure  based  on  theoretical  and  empirical  acoustical  models  for  the  estimation  of  the 
installation  effects  on  the  radiated  far  field  flyover  noise  emitted  from  a full  scale  jet 
engine.  Thus,  the  program  can  be  used  to  evaluate  the  total  noise  reduction  due  to  the 
installation  of  a noise  component  source  (inlet  fan,  exit  fan,  core,  turbine,  or  jet)  in  the 
proximity  of  a noise  shielding  surface. 

The  principle  considerations  of  storage,  input  and  output'  incorporated  in  the  design  of 
the  original  noise  estimation  program  developed  under  contract  NAS2-6969  have  been 
affected  by  the  new  additional  computer  software  developed  for  the  present  contract  as 
follows. 

1.  The  program  structure  ha3  remained  a two  level  overlay  with  the  primary  level 
controlling  the  reading  in  of  all  the  input,  printout  of  all  the  output,  and  linkage  to 
the  secondary  level.  The  secondary  level  includes  in  addition  to  the  subroutines  for 
flight  geometry  calculations,  extrapolation  corrections,  and  noise  component 
prediction,  new  packages  for  computing  the  sound  attenuation  due  to  wing 
shielding  of  jet  and  turbomachinery  noise  sources.  The  resultant  increase  in 
program  size  from  NAS2-6969  is  approximately  15%. 

2.  The  NAMELIST  statement  is  also  used  for  the  additional  input  data;  i.e.,  the 
engine  wing  geometry  data,  the  noise  component  shielding  parameters,  jet  edge 
interaction  data  and  jet  noise  shielding  data.  The  characteristics  of  NAMELIST 
provide  the  convenience  of  updating  the  data  set  from  case  to  following  case  by 
simply  specifying  only  the  single  or  array  variables  that  need  be  changed. 

3.  Output  is  printed  per  case  (a  case  is  all  data  for  a maximum  of  three  engine 
configurations  per  aircraft)  with  options  on  the  type  of  output  information  desired. 
The  format  and  details  of  the  output  are  identical  to  that  of  NAS2-6969  with  the 
addition  of  engine/wing  geometry  data  per  configuration,  wing  shielding  data  and 
predicted  corrections  for  shielding  in  each  noise  component. 

The  general  procedural  steps  for  processing  data  in  the  estimation  of  noise  including 
installation  effects  are  largely  a function  of  the  order  of  noise  component  inputs  which 
are  dependent  on  the  engine  (source)  configurations  for  a set  of  aircraft  conditions  (i.e., 
a case). 

These  aircraft  conditions  together  with  information  pertinent  to  the  current  case 
(described  under  general  data  parameters)  are  the  first  set  of  data  input  to  the  program 
for  a case.  The  primary  of  executive  overlay  initiates  the  flight  path/observer  geometry 
and  extrapolation  corrections  overlay,  to  calculate  the  aircraft  coordinates,  the  sound 
propagation  distance  and,  the  sound  transmission  time  for  each  observer  position  for  the 
17  angles  (10°,  20°,  . . . 170°)  between  the  flight  path  and  a line  to  the  observer. 
Extrapolation  corrections  including  sound  attenuation  due  to  spherical  divergence, 
atmospheric  absorp>>on,  extra-ground  attenuation,  and  ground  reflection  are 


accumulated  in  an  array  (array  size  = number  of  spectra  bands  times  17  aircraft 
positions  times  the  number  of  observer  positions)  used  for  extrapolating  the  index 
spectra  of  each  noise  component  and  the  total  index  spectra. 

If  shielding  is  indicated  in  the  general  data  parameters  for  the  particular  engine 
configuration,  then  the  engine/wing  geometry  data  for  that  configuration  must  be  input 
for  use  in  the  engine/edge  geometry  of  the  wing  shielding  package. 

From  hereon,  until  the  next  case,  there  could  be  up  to  three  engine  configurations 
where  the  inputs  for  each  configuration  relate  to  the  particular  noise  components 
making  up  the  configuration.  Data  for  each  configuration  consist  of  a title  card  foiled  by 
data  for  each  selected  noise  component,  followed  by  wing  shielding  data  for  the  noise 
component  if  shielding  is  indicated  for  the  particular  noise  component.  If  more  than  one 
case  is  to  be  input,  then  it  is  mandatory  that  the  order  and  number  of  noise  components 
in  each  engine  configuration  be  maintained  for  all  cases.  After  each  noise  component 
data  and  associated  shielding  data  is  read  in,  the  executive  module  links  the 
corresponding  overlay. 

For  engine-over-wing  (EOW)  configured  aircraft,  the  noise  source  components  of 
relevance  are  primary  and  secondary  jet,  core,  turbine,  and  inlet  and  exit  fan. 
Considering  these  noise  sources,  computation  modules  are  provided  to  execute  with 
some  variations  the  following  sequential  computation  steps. 

a.  Calculation  of  the  directivity  angles,  the  projected  directivity  angles  and  the 
projected  elevation  angles  for  each  sideline  observer  position  and  noise  component 

b.  Calculate  lining  attenuation  spectra  if  applicable 

c.  Calculate  noise  prediction  spectra  without  effects  of  shielding 

d.  Compute  free-field  index  spectra 

e.  Include  multiengine  effects 

f.  Convert  spectra  to  full  octave,  if  applicable 

g.  Calculate  attenuation  due  to  shielding 

h.  Sum  configuration  effects,  attenuation  due  to  shielding,  and  subtract  the  results 
from  the  bare-engine  index  spectra 

i.  Compute  total  free-field  noise  spectra 

j.  Subtract  extrapolation  corrections  from  free-field  index  spectra 

k.  Compute  human  response  measures;  i.e.,  perceived  noise  level  (PNL),  tone 
corrected  PNL,  and  effective  PNL 


3 


to 


1.  Save  the  index  and  observed  sound  spectra  on  scratch  data  files 


f 


After  all  noise  components  for  all  the  configurations  in  a case  have  been  processed,  the 
total  free-field  index  spectra  are  extrapolated  and  the  human  response  measures  are 
calculated  based  on  the  total  extrapolated  spectra.  Depending  on  input  specifications  for 
the  output  reports  desired,  information  on  the  scratch  data  files  is  printed.  The  program 
will  continue  to  process  input  data  until  there  are  no  more  cases;  i.e.,  the  program 
terminates  execution  whenever  an  end-of-record,  end-of-file,  or  end-of-information  code 
is  encountered  on  the  input  file. 
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3.0  INTERFACE  WITH  NAS2-6969  CONTRACT  SOFTWARE 


The  interface  of  the  original  NAS2-6969  contract  computer  software  with  the  new 
software  developed  for  contract  DOT-FA72WA-3497  has  resulted  in  minor  modifications 
of  the  original  software.  These  modifications  are  mainly,  but  not  exclusively  limited  to 
the  noise  source  component  modules  effected  by  shielding  and,  the  executive  routine. 
The  interface  is  structured  so  that  all  references  to  wing  shielding  can  be  completely 
excluded,  thus  retaining  the  external  operational  characteristics  of  the  NAS2-6969 
contract  software.1  Figures  A-l  and  A-2  (similar  to  fig.  1 and  3 of  NASA-CR-1 14650) 
illustrate  the  overall  structure  of  the  noise  prediction  program  including  interface, 
additional  input  and  output,  and  noise  source  component  modules  influenced  by  wing 
shielding.  Figure  A3  is  a macro  flow  chart  of  the  subroutine  SHLDSP,  which  interfaces 
NAS2-6969  software  with  the  evaluation  of  turbomachinery  noise  shielding  effects.  That 
is,  noise  source  prediction  modules  INLET  FAN,  EXIT  FAN,  CORE,  and  TURBINE  are 
linked  with  the  wing  shielding  module  SHATTN,  through  the  routine  SHLDSP  which 
returns  shielded  SPL  spectra  at  flight  index  conditions  for  each  of  <le  turbomachinery 
noise  component  estimation  modules.  The  total  programming  add  ion  to  NAS2-6969 
generates  an  increase  in  program  field  length  of  approximately  1 5'v •. 

Jet  noise  shielding  modules  for  the  primary  and  primary  plus  secondary  jet  noise 
sources  are  made  up  of  the  subroutines  JNSA  and  JNSASP.  The  shielding  noise 
attenuation  is  returned  to  the  jet  noise  source  prediction  module,  JET  (fig.  A-4). 

3.1  WING  SHIELDING  PACKAGE* 

The  wing  shielding  package  may  be  considered  a system  of  computer  programs  designed 
to  simulate  a method  of  quantitatively  estimating  the  flyover  noise  reductions 
attainable  from  engine  placement  in  the  proximity  of  wing  structures  which  provide 
shielding.  The  package  of  computer  programs  to  predict  shielded  noise  is  common  to 
engine  component  noise  sources  fan  (inlet  and  exit),  core,  and  turbine,  (see  fig.  A-2  and 
A-5).  The  jet  noise  source  component  (primary  and  secondary)  is  treated  separately  with 
its  own  special  shielding  programs  (see  fig.  A-2  and  A-6). 

3.2  CALCULATION  SEQUENCE 

Figure  A-5  is  an  outline  of  the  subtasks  that  contain  the  sequence  of  calculations  to 
determine  wing  shielding  attenuation.  The  module  SHATTN  and  its  externals  are 
common  to  noise  source  components  turbine,  core,  inlet  and  exit  fan.  The  module 
linking  SHATTN  to  the  noise  components  of  NAS2-6969  is  the  interface  routine 
SHLDSP  as  shown  in  figure  A-3.  The  calculation  sequence  for  determining  jet  noise 
shielding  is  shown  in  figure  A-6;  the  linkage  to  the  noise  component  through  JNSA  is 
shown  in  figure  A-4. 


’New  additions  to  contract  NAS2-69S9  software  due  to  KAA  contract  DOT-FA74WA-3497  such  as 
subroutines,  data  or  data  sets  are  indicated  with  an  * , and  software  with  onl\  pertinent 
modifications  is  indicated  with  an  **. 
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Figure  A ■ 1. -Macro  Flow  Chart  of  The  Executive  Routine- 
Showing  Integration  of  Noise  Components 
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calc:  Corrections  for  spherical  divergence, 
atmospheric  absorption,  extra  ground  attenua- 
tion, ground  reflection 


I Figure  A-1.~(Con  tinued) 

, 
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Test  noise  component  type  and  branch  to  specified  noise 
component. 


ITY  PE= 1 
Primary , jet** 


ITYPE=3 

Core  & turbine** 


ITYPE=4 
Inlet  fan** 


ITYPE=5 

Exit (aft) fan** 


ITYPE=6 

Augmentor  wing 


ITYPE=7 
Blown  flap 


ITYPE=8 
Lift  fan 


ITYPE=9 

Ejector  suppressor 


ITYPE=1Q 

Propeller 


ITYPE=1 1 
Helicopter 


ITYPE=12 
Measured  data 


ITYPE=1 3 
Jet  edge  * 


Transfer  to  calculate 


next  noise  componen 


PNLSUB 

calc:  Perceived  noise  level,  effective  perceived 
noise  level,  tone  corrections  from  total 
predicted  noise 


NO I SO** 

print:  observed  spectra  for  total  noise  of  all 

components  and  all  report  types  requested 


Figure  A- 1.  —( Concluded ) 
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FLa  =*  OB  to  105136B 


OVERLAY  0 (Main  Overlay  Program  TEE33QA**) 

Purpose: 

Routines 

called: 

0 

Initialization 

ABSORP 

ANGLES 

CONVR 

CORSPL** 

EPNLI 

k ERROR 

EWGIO* 

EWGOUT* 

FAAISO 

INITL 

0 

Input/Output 

L INPUTG** 

LININS 

N INPUT** 

NOISO** 

NSRIO* 

PND64 

PNLSUB 

PRINTH 

PRTSH1* 

SN 

0 

S.I.  or  English 

SEARCH 

TBLU1 

TBLU2 

TERPI 

TERP2 

unit  conversion 

TONCOR 

T0NES3 

URPNL 

WSOSTO* 

XPM1DX 

0 

Calls  Overlays 
1 , 2,  ...  or  6 

f 
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FLa  * 105137B  to  112756B 


OVERLAY  1 (Secondary  Overlay  Program  ONE) 

Purpose: 

Routines 

called: 

o Subroutine  FLTGEO  solves 

flight  geometry  and  determines 
atmospheric  properties. 

ATMOSP 

ERROR 

SORTX 

TBLUI 

o Subroutine  NEXTCR  calculates 
spherical  divergence  attenua- 
tion, atmospheric  absorption, 
extra  ground  attenuation,  and 
ground  reflection  interference 
corrections. 

AVGALF 

ATMOSP 

TBLUI 

ABSORP 

CPOLAR 

EGACAL 

DPSD 

GRDRFX 

FW 

FLa  « 1Q5137B  to  141024B 


OVERLAY  2 (Secondary  Overlay  Program  TWO**) 

Purpose: 

Routines 

calledb: 

0 

Subroutine  AFT**  calculates  '1 

BUZSAW 

CFI* 

CPOLAR 

EDGEDI*- 

exit  fan  noise  component 

EEGEOM* 

EOGEOM* 

ESHLDG 

FANNOS 

FANPED 

REFRACT* 

RESCAL 

SDSPLP* 

0 

Subroutine  INLET**  calculates 

SHATTN* 

SHELL 

SHELXH 

SHLDSP* 

Inlet  noise  component 

UNIFLW* 

UNIT 

VECN* 

DOTP*  ' 

J 

CROSP* 

WSHOUT* 

ZERO 

0 

Subroutine  UET**  calculates 

ESHLDG 

JETNOS 

JETPED 

JNSA* 

jet  noise  component 

JNSASP* 

TBLU2 

TBLU3 

TERP2 

TERP3 

UNIT 

ZERO 

a)  Field  length  (FL)  noted  Includes  system  routines. 

b)  All  noise  source  component  modules  also  call  the  following  routines 
TrTaddltion  to  those  noted  In  the  blocks;  ANGLES,  CORSPL,  LGMTRY, 
LINCOR,  LINING,  LININS,  NOISO,  PNLSUB,  PWRSUM,  SEARCH,  SORTX,  TBLU1, 
and  TERP1 . 


r 
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FLa  = I05137B  to  124074B 


OVERLAY  3 (Secondary  Overlay  Program  THREE**) 


Purpose: 

Routines 

called6: 

0 

Subroutine  COREN**  calculates 

CFI* 

CPOLAR 

EDGED1* 

EEGEOM* 

core  and  turbine  noise 

EOGEOM* 

ESHLDG 

REFRACT* 

SDSPLP* 

components . 

SHATTN* 

SHELL 

SHELXH 

SHLDSP* 

UNIFLW* 

WSHOUT* 

VECN* 

DOTP* 

CROSP* 

0 

Subroutine  EJECT  calculates 

ESHLDG 

MENOZZ 

TBLU2 

TERP2 

jet  noise  component  for  an 
ejector/multi -element  suppres- 

ZERO 

so  r. 

0 

Subroutine  SPECAN  calculates 
jet  noise  component  for  an 
augmentor  wing/slot  nozzle. 

TBLU2 

TERP2 

ELa  ° 105137B  to  115707B 


OVERLAY  4 (Secondary  Overlay  Program  FOUR) 

Purpose: 

Routines 

called6: 

o Subroutine  COPTR  calculates 
helicopter  or  tilt  rotor 
noise  component. 

BESJ 

ERROR 

JBES 

o Subroutine  LIFTFN  calculates 
lift  fan  noise  component. 

BUZSAW 

RESCAL 

ESHLDG 

UNIT 

FANNOS  FANPED 

ZERO 

o Subroutine  PROP  calculates 
propeller  noise  component. 

ERROR 

FLa  * 105137B  to  135123B 


OVERLAY  5 (Secondary  Overlay  Program  FIVE) 

Purpose: 

Routines  called6: 

o Subroutine  MEASRD  interpolates 

TBLU3  TERP3  TERP2 

on  tabular  Input  of  measured 

noise  data. 

. FLa  =■  105137B  to  113643B 

OVERLAY  6 (Secondary  Overlay  Program  SIX*) 


Purpose:  Routines  called6: 

o Subroutine  BLWFLP  calculates  ESHLDG 

externally  blown  flap  noise  TERP2 

component. 

o Subroutine  JEINT*  calculates  ZERO 

jet/edge  Interaction  noise  for 
engine-over-wing  configurations. 

TBLU3 

TBLU2 

TERP3 

Figure  A -2.-( Concluded) 
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Figure  A-3.— Routine  SHLDSP  to  Interface  NAS2-6969*  With  Wing 
Shielding  Calculation  Module  (SHA  TTN) 


II 


inni^ 


ENGINE/EDGE  GEOMETRY 


OBSERVER  SIDELINE  DISTANCE 


EDGE/OBSERVER  GEOMETRY 


JET  FLOW  REFRACTION  EFFECTS 


UNIFORM  FLOW  EFFECTS 


EDGE  DIFFRACTION 


SUPERPOSITION  OF  EDGE  SOLUTIONS 


COMPUTE  SHIELDING  ATTENUATION 


Figure  A-5.-Calculation  Sequence  in  Wing  Shielding  Module  SHATTN * 
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3.2.1  EXIT  FAN,  CORE,  AND  TURBINE  NOISE  COMPONENT 


Figure  A-7  shows  the  general  logic  flow  of  the  calculations  for  the  noise  components; 
exit  fan  (AFT),  turbine  and  core  (COREN).  When  the  effect  of  shielding  is  to  be 
included,  noise  components  core  and  turbine  must  be  run  separately  due  to  the  unique 
noise  component  wing  shielding  data  for  each  component. 

3.2.2  INLET  FAN  OR  COMPRESSOR  NOISE  COMPONENT 

Figure  A-7  shows  the  general  logic  of  the  calculations  for  the  noise  component  inlet  fan 
or  compressor  (INLET).  The  programming  steps  are  similar  to  the  discharge 
turbomachinery  noise  components,  except  when  shielding  is  included  there  are  no  jet 
flow  refraction  effects,  thus  reducing  the  input  requirements  for  wing  shielding  data  for 
the  inlet  fan  noise  component. 

3.2.3  JET  NOISE  COMPONENT 

The  procedure  to  predict  the  jet  noise  shielding  attenuations  for  an  EOW  configuration 
is  unique  to  the  noise  components  primary  and  primary  plus  Stecondary  jet,  and  is 
therefore,  treated  separately  from  noise  components  inlet  fan,  aft  fan,  core, and  turbine. 
The  attenuations  are  to  be  subtracted  from  the  predicted  unshielded  jet  noise  spectra 
assumed  to  be  at  index  conditions.  Figure  A-4,  given  previously,  shows  the  logic  for 
subroutine  JET  with  the  additional  subroutine  (JNSA)  necessary  for  jet  noise  shielding. 

3.2.4  JET/EDGE  INTERACTION  NOISE  COMPONENT* 

The  jet/edge  interaction  noise  component  is  a new  noise  source  that  must  be  considered 
in  evaluating  an  EOW  aircraft  configuration.  This  noise  is  caused  by  the  interaction  of 
the  jet  exhaust  with  the  wing’s  trailing  edge.  A computer  module  has  been  added  to  the 
NAS2-6969  contract  software  for  estimation  of  this  type  (ITYPE  = 13)  of  noise. 
Figure  A-8  shows  the  computation  steps  for  calculations  done  by  this  new  noise  source 
estimation  module. 

The  jet/edge  interaction  noise  source  module  is  linked  to  the  main  program  in  a manner 
similar  to  that  used  for  the  other  noise  source  estimation  modules  available  from  the 
NAS2-6969  contract.  Therefore,  the  input  scheme  for  this  noise  component  uses  the 
same  NAMELIST  name;  i.e.,  $NOISIN,  as  used  for  the  data  sets  of  the  other  noise 
components. 


Figure  A-7.  —General  Logic  Flow  of  Noise  Components— Aft  Fan,  Core, 
and  Turbine,  and  Inlet  Fan 
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4.0  INPUT  DATA 


I 


t 
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In  order  to  facilitate  the  use  of  this  program,  many  of  the  data  inputs  have  prestored 
default  values  in  the  program  and  are  defined  in  the  appropriate  parameter  description 
section.  An  examination  of  the  following  sections  will  indicate  the  necessary  input  data 
needed  for  the  various  options  built  into  the  program;  e.g.,  establishing  different 
atmospheric  conditions;  using  calculated  ground  reflection  corrections  or  a 3-dB  default 
value;  using  measured  data  to  define  a noise  configuration  instead  of  that  obtained  by  a 
prediction  method;  specifying  whether  shielding  is  to  be  included  for  a particular 
configuration.  The  data  input  procedure  is  fully  described  in  the  sections  following  and 
includes  in  addition  to  data  for  shielding  all  the  input  as  documented  in  NASA 
CR-1 14650.  Table  A-l  presents  a cross  reference  of  the  engineering  and  programming 
documents  by  noise  type  and  section  number. 

The  program  processes  one  or  more  cases  of  data  in  a single  run.  A case  consists  of  from 
one  to  three  different  types  of  propulsion  systems  or  configurations.  The  data  for  each 
configuration  are  described  in  terms  of  a subset  of  noise  components.  The  order  of  input 
for  each  case  is  fixed  after  the  first  case.  Inputs  which  remain  constant  throughout  the 
run,  need  no^  be  redefined  for  succeeding  cases. 

The  first  input  per  case  consists  of  an  80  column  alphanumeric  title  record  which  will 
be  output  as  part  of  the  heading  describing  the  case.  If  the  same  heading  is  desired 
throughout  the  run,  only  blank  records  need  be  entered  in  succeeding  cases.  The  title 
record  is  followed  by  data  cards  describing  the  general  conditions  which  apply 
throughout  the  case.  These  cards  are  coded  in  a FORTRAN  NAMELIST  format.  The 
NAMELIST  format  requires  that  the  first  records  start  in  column  2 with  an  initiator 
name  SGDATA  on  General  Data  Parameters.  The  order  of  the  parameters  is 
independent  but  must  have  the  variable  name,  followed  by  an  equals  mark,  followed  by  >. 

the  assigned  value.  The  parameters  must  be  separated  by  commas;  blanks  are  ignored. 

Any  succeeding  cards  must  start  in  or  after  column  2 and  may  go  though  to  column  80. 

A terminator  $ must  follow  the  iast  parameter  in  each  data  set. 

A variable  name  can  be  followed  by  a single  value,  or  a series  of  values  in  the  case  of 
variables  which  require  an  array  of  data.  The  general  input  format  is  as  follows: 

Variables  Having  Single  Values 
Example:  integer  type  decimal  type 

NENG  = 3 SLOPE  = .25 

Variables  Which  Have  Multiple  Values 
Example:  Filling  the  complete  array  for  a 10  element 

array 

SLDIST(l)  or  SLDIST  = 100.,  200.,  150..  175..  500.,  750.,  1000., 

1200..  1500.,  2000., 

Filling  part  of  the  array 
SLDIST03)  - 300.,  400., 
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Each  SN0ISIN  record  specifies  an  ITYPE  module  and  all  arguments  in  the  record  are  for  that 
module  (lining  and  configuration  options  are  included  in  the  record  as  required). 


Example  of  NAMELIST  GDATA: 

$GDATA  ALTOG=  100.,  ALTPG  = 20t)..  SLOPE=.l.  AMACH  = .15.  NOBS=2, 
SLDIST(l!=100„  500.,  ISPTRM=  1 . NTENG  = 2,  INSE0W<1)-1.  1,1$ 
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Sample  cases  with  input  data  are  presented  in  this  appendix  A,  sections  5.1  and  5.2. 

If  shielding  is  specified  in  $GDATA,  then  engine  wing  geometry  data  SEWDATA  must 
be  input  as  the  next  data  set  for  the  particular  configuration.  As  the  geometrical 
relationship  for  the  wing  and  each  propulsion  system  may  vary,  the  $EWI)ATA  data  set 
is  input  for  each  configuration. 

Following  the  SGDATA  data  set  or  when  shielding  is  specified  the  $EWDATA  data  set 
is  an  80  column  alphanumeric  title  record  describing  the  first  noise  source 
configuration.  Separate  data  sets  using  the  NAMELIST  initiator  $NOISIN  describe  each 
noise  component  in  a configuration.  If  more  than  one  configuration  is  to  be  included,  the 
preceeding  type  of  input  is  included  for  each  additional  configuration  starting  with  the 
title  record  and  ending  with  the  SNAMELIST  terminator  for  the  last  noise  component 
in  the  $N()ISIN  data  set  As  mentioned  previously,  there  is  a maximum  of  three 
different  types  of  powerplants  per  run;  therefore  a case  may  be  processed  having 
different  engine  configurations  using  any  or  all  of  the  13  different  noise  modules 
depending  on  the  nature  of  the  particular  configuration.  Reference  to  the  engineering 
analysis  is  itemized  in  the  engineering  programming  cross  reference  shown  in  table  A-l. 
Each  noise  component  data  set  is  described  in  the  following  sections.  A special  input 
variable  terminator  is  not  required  to  end  the  run  since  an  end-of-record  check  is 
repeatedly  made  on  the  input  data  fiie  at  the  end  of  each  case 

If  shielding  is  considered  and  the  noise  component  is  effected  bv  shielding,  the  i the 
SNOISIN  data  set  is  followed  bv  $SHLI)AT  data  set  describing  the  shielding  data 

When  using  the  program  tor  multiple  muse  source  configurations,  an  input  value  for 
NTENG  2 or  is  required  in  the  $GI>AIA  set  When  more  tnan  one  type  of  noise 
component  is  used  in  a configuration  an  input  NTYPF  m where  m is  the  total  noise 
components  for  the  configuration  i*  needed  in  tin  first  $N<>1S1N  data  set  for  that 
configuration  The  setem  ...wring  the  sample  a*.  *eiti.-i,  -h-mld  assist  in  clarifying 
(he  tnpu*  data 
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defined  by  input  or  by  the  default  value(s).  In  order  to  facilitate  finding  the  description 
of  a parameter,  the  list(s)  of  inputs  within  each  of  the  following  sections  have  been 
alphabetized  with  respect  to  the  names  of  the  variables. 


Description 


Airport  altitude. 

Aircraft  height  above  ground  at  Y = 0 
(fig.  6 of  NASA  CR-1 14649). 

Observer  height  above  the  ground. 
Aircraft  Mach  number.. 

Number  of  perceived  noise  decibels 
down  from  maximum,  used  to  determine 
the  integration  interval  for  the 
EPNL  calculations. 

Pressure  of  homogeneous  atmosphere 
defined  by  user.  (Note:  Used  only 
if  I ATMOS  = 4.) 

Relative  humidity  of  homogeneous 
atmosphere  defined  by  user.  (Note: 

Used  only  if  IATMOS  =4.) 

Temperature  of  homogeneous  atmosphere 
defined  by  user.  (Note:  Used  only 
if  IATMOS  = 4.) 

Constant  percent  relative  humidity  delta 
that  is  added  to  ISA.  (Note:  Used 
only  if  IATMOS  =1.) 

Constant  pressure  delta  that  is 
added  to  ISA.  (Note:  Used  only  if 
IATMOS  = 1.) 

Constant  temperature  delta  that  is 
added  to  ISA.  (Note:  used  only 
if  IATMOS  = 1.) 

Engine  performance  parameter  used 
to  correlate  data  output  for  the 
Noise  Contour  Program  or  to  specify 
operating  condition  when  using  the 
"measured  data"  module. 

Data  frequencies  for  the  ground 
impedance  data  curve  used  in  the 
ground  reflection  calculation.  The 
number  of  values  to  be  specified  is 
ND.  (Note:  Used  only  if  IGDR  = 0.) 


Variable 

Name 

Units 

Default 

AALT 

M(ft) 

0. 

ALTPG 

m(ft) 

0. 

ALTOG 

m(ft) 

0. 

AMACH 

BCG 

PNdB 

10. 

CPRES 

Atm(psia) 

CRHUMD 

%RH 

CTEMP 

K(°R) 

DHUMID 

%RH 

0. 

DPRES 

Atm(psia) 

0. 

DTEMP 

K(°R) 

0. 

EPP 

0. 

FLD(l) 

Hz 

FLD(25) 

Hz 
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Variable 

Name 

Units 

Default 

Description 

FLR 

PNdB 

.90. 

Noise  floor  for  EPNL  calculations.  , 

Values  less  than  FLR  are  not  included  ' 

when  computing  the  duration  for  EPNL. 

IAIR 

0 

Specifies  whether  the  air  absorption 
coefficients  are  calculated  by  the 
program  or  defined  by  the  user.  Set 
equal  to: 

0 if  program  calculates 

1 if  user  defines  (see  UAIRAB) 

-1  if  program  retains  coefficients  \ 

from  previous  case  ’ 

IATMOS 

0 

Specifies  the  type  of  atmospheric 
conditions  used  by  program.  Set 

equal  to:  j 

0 for  nonhomogeneous  international 
standard  atmosphere  (ISA) 

1 for  nonhomogeneous  ISA  plus  user  ' 

defined  constants  for  relative 

IDOP 


humidity,  pressure,  and 
temperature  added  to  ISA.  (Note: 

See  DHUMID.  L)PRES.  DTEMP.) 

2 for  nonhomogeneous  atmospheric 
conditions  that  are  defined  bv 
the  user.  (Note.  See  NTEMP, 

TALT.  TEMP.  NPRES,  PALT,  PRES. 
NHUMIl).  RALT.  RHUMID.) 

3 homogeneous  atmosphere  of  1 
ATM  = 14.696  psia;  288.16  K 
(15°  C>  - 518.688"  R (59"  Ft; 
im  RH 

4 homogeneous  atmosphere  defined 
by  user.  (Note  See  CPRES, 
CRHUMD,  CTEMP  ) 

Doppler  shift  switch. 

0 for  no  Doppler  shift  factor 

1 Doppler  shift  factor  included 
(Note:  For  fan  noise  the 
correction  is  for  frequency 
only.) 

2 Doppler  shift  factor  included 
(Note:  For  fan  noise  a correction 
is  made  for  both  the  frequency 
and  the  noise  levels.) 
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Variable 

Name 


Units 


Default 


Description 


i ' 

1 1 
i 


I EGA 


0 


IGDR 


1 


*INSEOW(  1 ) 0 

INSEOW<3>  0 


IOUT(l)  0 

IOUT(7)  0 


(Note:  App.  A,  sec.  5.0  output 
data  illustrates  the  various  output 
reports  described  by  this  input.) 


Specifies  whether  corrections  for 
extra  ground  attenuation  are  to  be 
applied  while  extrapolating  the 
noise  level  from  the  airplane  to 
the  observer.  Set  equal  to: 

0 if  EGA  desired 

1 if  EGA  is  not  wanted 

Specifies  whether  corrections  for 
ground  reflection  based  on  +3-dB 
delta,  or  calculated  corrections 
(app.  A.  NASA  CRl  14649*  are  to  be 
applied  while  extrapolating  the  noise 
level  from  the  airplane  to  the 
observer.  Set  equal  to: 

0 if  calculated  corrections  are 
made  (Note:  See  XKN.  ND.  ELI). 

ZNR.  ZN1) 

1 if  a delta  of  +3  dB  is  to  be 
used 

An  array  specifying  whether  shielding 
is  to  be  included  for  a particular 
engine  configuration  in  which  case 
engine  wing  geometry  data  i$EWl)ATA 
data  set  i is  to  be  input  for  that 
configuration. 

0 engine  wing  geometry  data  are 
not  to  be  input  (unshielded* 

1 engine  wing  geometry  data  must 
be  input  (shielded) 

An  array  to  indicate  the  output  reports 
desired.  The  selection  can  include  seven 
different  output  reports  or  a default 
report  which  gives  a heading.  PNL.  TCPNL. 
time  array.  EPNL,  and  a one  page  summary  of 
assumptions  under  which  the  run  was  made. 
Order  of  input  is  immaterial  e g.. 

IOUTi  1 ) - 5,  4.  3,  2.  1 functions  the  same  as 
IOUT(l)  =1,2.  3.  4.  5 
= 1 Type  1 report;  total  SPL  at  the 
observer  for  each  of  the  17  angles 

(10°.  20".  30" 170")  for  each  of 

the  frequency  bands  and  observer  positions 
= 2 Type  2 report;  summary  of  options  assump- 
tions under  which  the  case  was  made. 


■ 1 

•.  I 
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Variable 

Name 


Units 


Default  Description 


u 


b 

I 


M 


j 

I 


H 

r 


i ' 

[ i 
f i 

t. 

i . 

V 

' 


*ISFE 


ISPTRM 


IUNIT 


ND 


NHUMID 


= 3 Type  3 report;  SPL.  for  each  component, 
at  the  observer  in  the  same  manner  as 
type  1 . 

= 4 Type  4 report;  flightpath  observer  geom- 
etry and  engine  wing  geometry  if 
applicable, 

= 5 Type  5 report;  extrapolation  corrections. 
= 6 Type  6 report;  total  free-field  index 
(radius  of  1 mi  SPL's  for  all  angles, 
frequencies  and  observer  positions. 

* 7 Type  7 report;  free-field  index  spectra 
for  each  component.  If  shielding  is 
specified,  the  type  7 report  selection 
produces  output  of  wing  shielding  com- 
ponent data. 

0 Specifies  if  flight  effects  per 

uniform  flow  theory  are  to  be  included 
in  wing  shielding  calculations. 

1SFE  = 0 for  inclusion  of  flight 
effects 

ISFE  * 0 for  prediction  based  on 
static  test  results 

0 Specifies  the  type  of  frequency  bands 

to  be  used  in  the  calculations.  Set 
equal  to: 

0 for  24  preferred  1 3 octave  bands 

1 for  8 preferred  1/1  octave  bands 

0 Specifies  whether  input  parameters  and 

output  reports  are  in  S.l.  or  English 
units. 

0 = S.l. 

1 = English  units 

3 Number  of  data  points  for  the  ground 

normalized  complex  impedance  curve. 

(3  « ND  =s  25)  (Note:  Use  only  if 
IGDR  = 0.  see  XKN,  FLD,  ZNR.  ZNI.) 

Specifies  the  number  of  entries  in  each 
of  the  percent  relative  humidity  (RHUM1I) 
versus  RALT)  tables  that  are  defined  by 
the  user  for  nonhomogeneous  atmospheric 
conditions.  (2  --  NHUMID  **  50)  (Note: 

Use  only  if  I ATMOS  = 2.) 


24 
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Variable  Units  Default  Description 

Name 


NLOPT  0 


NOBS  1 

NPRES 


NTENG  1 


NTEMP 


PALT(l)  mi  ft  i 

PALTI50)  ml  ft) 

PRES(l)  Atmlpsial 

PRES(50)  Atm(psia) 


Specifies  table  output  for  noise 
contour  estimation  on  file  TAPE20. 

0 - no  output 

! - F PNL  versus  KPP.  elevation 

1 ngle.  log , o of  the  range  at  CPA 
2 - same  except  peak  PNL 
(Note:  Output  reporting  must  be  set 
to  TYPE  1 or  default-see  lOl'T  If 
NLOPT  0.  each  case  will  have  a noise 
level,  an  engine  performance 
parameter,  an  elevation  angle,  and 
the  log | <(  of  the  off-axis  range 
written  to  file  TAPK2 0.) 

Number  of  observe-,  defined  in 
SLD1ST  table  i ! NOBS-  10.) 

Specifies  the  number  of  entries  in 
each  of  the  pressure  iPRES>  versus 
altitude  iPALT>  tables  that  the  user 
defines  for  nonhomogeneous  atmospheric 
conditions.  <2  * NPRES  ’ ”>0>  (Note: 

Used  if  1 ATMOS  2.i 

Specifies  the  number  of  distinct 
tvpes  of  engine  configurations  to 
be  considered  (NTENG  * 3).  Noise 
component  parameters  must  be  defined 
for  each  different  source. 

Specifies  the  number  of  entries  in  each 
of  the  temperature  (TEMP)  versus 
altitude  (TALTi  tables  that  have  been 
defined  by  the  user  for  nonhomogeneous 
atmospheric  conditions.  (2  - NTEMP  5(11 
(Note:  Used  only  if  IATMOS  - 2.) 

Each  entry  in  this  table  defines  the 
altitude  for  the  pressure  defined  by 
the  corresponding  entry  in  the  PRES 
table.  (Note:  Used  only  if 
IATMOS  = 2.  see  NPRES.) 

Each  entry  in  this  table  defines  the 
pressure  for  the  altitude  defined  by 
the  corresponding  entry  in  the  PALT 
table.  (Note:  Used  only  if 
IATMOS  - 2.  see  NPRES.) 


25 


I 


1 


I 


I 

1 


I 


Variable 

Nairn* 

Units  Default 

Description 

RALT(l) 

mtft) 

Each  entry  in  this  table  defines  the 
altitude  for  the  relative  humidity 

RALT(50> 

mift ) 

defined  by  the  corresponding  entry 
in  the  RHUMID  table.  (Note:  Used 

only  if  IATMOS  = 2.  see  N HUMID. ) 

RHUMIDtl) 

QRH 

Each  entry  in  the  table  defines  the  percent 
relative  humidity  for  the  altitude 

RHUMIDioOt 

ci  RH 

defined  by  the  corresponding  entry  in 
the  RALT  table.  (Note:  Used  only  if 
IATMOS  = 2,  see  NHUMID.) 

SLDIST(l) 

m(ft> 

Sideline  position  of  1st  observer 
(see  NOBS). 

SLDIST(N) 

ml  ft) 

Sideline  position  of  Nth  observer. 

SLOPE 

0 

Aircraft  climb  gradient.  (Tangent  of 
climb  angle,  i 

TALTUi 

mtftt 

Each  entry  in  this  table  defines  the 
altitude  for  the  temperature  defined 

TALT(50) 

m(ft) 

by  the  corresponding  entry  in  the 
TEMP  table.  (Note:  Used  only  if 
IATMOS  2.  see  NTEMP.i 

TEMPO ) 

K(°R) 

Each  entry  in  this  table  defines  the 
temperature  for  the  altitude  defined 

TEMPI  50) 

K(°R) 

by  the  corresponding  entry  in  the 
TALT  table.  (Note:  Used  only  if 
IATMOS  -=  2.  see  NTEMP.i 

TCG 

sec  10. 

Normalizing  time  constant  in  seconds 
used  in  the  EPNL  calculations  (Note: 

See  BCG.  FLR. > 

UAlRABl  1 ) 

dB  kmldB  1000  ft) 

User  defined  air  absorption  coefficient 
for  the  frequency  bands.  'Note.  Used 

UAIRAB(N) 

dB  kmldB  1000  ft) 

only  if  I AIR  = 1).  N - 8 for  1 1 O.B.; 
N - 24  for  1 3 O.B. 

XKN 

Wave  number  ratio.  XKN  K K,,  or  (': 

where  C„  = speed  of  sound  in  air: 

C - speed  of  sound  in  ground. 

(Note:  Used  only  if  IGDR  0.  see  FLO. 
ND.  ZN1,  ZRN.i  RESTRICTION:  XKN  <>. 

ZNR(l) 

(See  app.  A, 

Real  part  of  (Z,  Z0)  for  normalized 

NASA  CR114649.) 

ground  impedance  data  curve.  (Note: 

ZNRI25) 

ZNR  • 0.  see  ELI).  ND.  XKN.  ZNT.  used  on 

if  IGDR  - O.i 
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Variable 

Name 


Units 


Default  Description 


ZNI(l)  Imaginary  part  of  -Zj/Z0  for  normalized 

ground  impedance  data  curve.  (Note: 

ZNI(25)  See  FLD,  ND.  XKN,  ZNR;  used  only  if 

IGDR  = 0.) 

(Note:  The  reactance  of  the  ground  is 
usually  capacitive,  hence  negative.  The 
option  here  permits  the  user  to  specify 
positive  values  which  are  treated  as 
capacitive  reactances.* 

4.2  ENGINE/WING  GEOMETRY*  ($EWDATA  DATA  SET) 

This  section  describes  the  engine  wing  geometry  data  to  be  input  for  each  shielded 
distinct  type  of  engine  configuration  (see  NTENG).  Successive  configurations  need  only 
reflect  the  change  in  the  '•ata  set.  See  figure  A-9  for  an  illustration  of  the  engine  wing 


geometry. 

Variable 

Name 

Units 

Default 

Description 

DDLD 

Nondimensional  engine  length,  L D 

DDSD 

Nondimensional  distance  between 
engine  centerline  and  wingtip.  S D. 

DDXOD 

0 

Nondimensional  distance  in  axial 
direction  from  nozzle  exit  plane 
to  point  on  top  of  w ing,  X,,  I). 

DDX1D 

Nondimensional  distance  in  axial 
direction  to  trailing  edge  from 
point  on  top  of  wing.  Xj  D. 

DDX2D 

Nondimensional  distance  in  axial 
direction  to  leading  edge  from 
point  on  top  of  wing.  X > I). 

DDYOD 

1 

Nondimensional  distance  normal  to 
engine  centerline  from  top  of  wing. 
Y„D. 

DDYlD 

Nondimensional  distance  normal  to 
engine  centerline  from  the  trailing 
edge  to  top  of  wing,  Yj  D. 

DDY2D 

Nondimensional  distance  normal  to 
engine  centerline  from  the  leading 
edge  to  top  of  wing.  Yj  D. 

Fiaure  A-9.- Engine /Wing  Geometry 


Variable 

Name 

Units 

Default 

Description 

DIANE 

mfft) 

0.3048 

Nozzle  diameter,  D, 

DIHED 

deg 

0. 

Wing  dihedral.  £*3. 

IES 

0 

Indicator  for  number  of  wings  to  be 
considered  in  shielding  calculations. 
IES  = 0 for  both  wings  on  near  and 
far  side  of  airplane 
IES  = 0 for  just  the  wing  on  the  nea 
side  of  airplane 

SWPLE 

deg 

0. 

Leading-edge  sweep  angle.  o_». 

SWPTE 

deg 

0. 

Trail ing-edge  sweep  angle, 

4.3  NOISE  COMPONENT  PARAMETERS  ($NOISIN  INPUT  DATA  SETS) 

As  mentioned  previously  the  program  estimates  noise  for  aircraft  equipped  with  one  of 
three  different  types  of  propulsion  systems.  Each  configuration  is  treated  us  a set  of 
noise  components,  each  of  which  will  have  a separate  SNOISIN  data  set  Each 
configuration  is  treated  independently  and  may  consider  any  subset  of  the  following  Id 
noise  component  modules: 

1.  Primary  jet 

2.  Primary  and  secondary  jet 

3.  Core  and  turbine.  •» 

4.  Compressor  and  fan  inlet,  tt 

5.  Fan  exit.  # 

6.  Augmenter  wing,  # 

7.  Blown  flap 

8.  Lift  fan. 

9.  Ejector  suppressor.  # 

10.  Propeller 

11.  Helicopter,  propeller,  and  tilt  rotor 

12.  Measured  data 

13.  Jet  edge  interaction* 

(Note:  denotes  lining  attenuation  option.) 


r 


\ 

i 


i ' 

* i 

■ i 


iiit  .1  of  these  components  has  meaning  only  in  respect  to  the  internal  order  in  the 
computer  program  and  a means  of  identifying  the  data  inputs  per  noise  component  (e.g., 
if  augmenter  wing  noise  is  to  he  predicted,  the  input,  ITYPE  = 6,  in  the  $N0IS1N  data 
set  indicates  to  the  program  that  this  data  set  applies  to  augmenter-wing  noise;  also 
TT6.  GAMA6.  DELT6.  etc.,  are  inputs  unique  to  the  data  set  when  ITYPE  - 6).  [For 
the  component  jet  edge  interaction  ITYPE  = 13.  all  data  are  input  through  an  array 
variable  name  EDGVAR.]X  Once  the  order  of  calling  the  noise  source  estimation 
modules  is  selected  for  the  first  case,  it  must  be  maintained  for  all  subsequent  cases. 

The  identifying  number  also  allows  the  lining  attenuation  and  configuration  correction 
inputs  to  be  described  in  a separate  section  which  is  referred  to  by  the  various  noise 
component  sections  rather  than  repeated  descriptions  in  each  section.  This  is  done  in 
the  following  manner  for  any  noise  component.  If  lining  or  configuration  corrections  are 
desired,  the  input  data  set  variable  names  are  appended  by  an  identifying  number  on 
the  end  of  the  variable  name  (e  g..  LIN’3  ---  1 icore  and  tubine  lining  attenuation!; 
LING  ~ 1 i augmenter  wing  lining  attenuation i;  IC’ORl  1 iprimarv  jet  configuration 
corrections!;  IGOR]  I 1 (helicopter  configuration  corrections)  etc.).  There  is  no  lining 
correction  option  for  component  jet  edge  interaction,  but  configuration  corrections,  if 
desired,  ( ICOR13  1 ) may  be  included. 

Three  data  inputs  for  each  distinct  source  noise  configuration  have  a special  place  in 
each  configuration  data  set  (i.e.,  they  must  be  defined  in  the  first  SN01S1N  data  set  in 
each  configuration!.  These  inputs  are; 
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t 


NTYPE  Specifies  the  total  number  of  noise  components  iNOISIN  data  sets)  in  a 
configuration  It  informs  the  program  to  accept  data  for  NTYPE  noise 
components.  3 

NENG  Specifies  the  number  of  identical  powerplants  on  the  aircraft.  > 

INSHLD  Specifies  wing  shielding  for  the  configuration.  The  computer  program 
expects  $SHLDAT  data  sets. 


4.3.1  LINING  ATTENUATION  AND  CONFIGURATION  CORRECTIONS 


This  section  describes  the  data  inputs  for  noise  components  t including  the  separate 
component  for  measured  data)  when  configuration  corrections  are  desired  in  the 
estimation  process  for  a particular  noise  component.  For  lining  attenuation  corrections 
this  section  applies  only  to  the  core  and  turbine,  compressor  and  inlet  fan.  exit  fan. 
augmenter  wing,  lift  fan,  and  ejector-suppressor  noise  modules.  The  input  variable 
names  described  below  differ  between  noise  components  only  by  the  appended  number 
described  in  the  introduction  of  the  Noise  Component  Parameters;  i.e  . ICORm.  where 
m ^ 1 . 2,  3,  4 . . . 1 2.  13,  for  the  particular  noise  component  i$N()lSlN  data  set). 
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Variable 

Name 


Units 


Default  Description 


ICORm  0 = 0 indicates  no  configuration 

corrections. 

=■  1 indicates  AdB  corrections  are  a function 
of  directivity  angle  only 
= 2 indicates  AdB  corrections  are  a function 
of  frequency  ' 1 3 or  1 1 octave'  and 
directivity  angle 

(Note:  Changing  ICORm  from  1 to  2 or 
vice  versa  is  not  permitted  for  the 
same  noise  component  of  an  engine 
configuration.  Each  engine  configur- 
ation can  be  treated  independently  of 
the  other  by  use  of  NTENG  option  il 
this  comparison  is  desired.) 


LINm 

0 

= 0 indicates  no  lining  attenuation 
in  corrections. 

= 1 indicates  lining  attenuation  corrections 
are  calculated  by  program. 

Lining  Attenuation  Parameters: 

LINm  / 0 

CFm 

m/s(fps) 

Speed  of  sound  in  the  flow. 

EDHm 

ELOHm 

FMm 

m(ft) 

Effective  duct  height  for  lining. 
(Note:  Used  only  if  LGMm  = 0.) 

Ratio  of  effective  lining  treatment 
length  to  duct  height.  (Note: 

Used  only  if  LGMm  - 0.' 

Mach  number  of  the  flow.  (Note: 
FMm  is  negative  for  inlet  lining.' 

IDPm 

2 

Lining  design  point  option. 

- 1 for  single  design  point 

- 2 for  multiple  design  point 

ILAYm 

1 

- 1 for  single  layer  lining 

- 2 foi  double  layer  lining 

IMAm 

0 

Specifies  whether  program  calculates 
or  user  defines  the  peak  attenuation 

for  each  target  frequency. 

0 program  calculates 
■ 1 user  defines  PLAm  values 

Specifies  whether  program  calculates 
peak  attenuation  using  lining  geometry 
or  user-defined  effective  duet  height 


LG  Mm 


() 


Variable 

Name 


Units 


Default  Description 


j. 

I: 


c 


1 


and  ratio  of  treatment  length  to 
effect  duct  height. 

- 0 user  inputs  EDHm  and  ELOHm 
= 1 user  inputs  lining  geometry 
(See:  RlWm  and  TLm.) 

(Note:  use  only  if  IMAm  = 0.) 

NTFm  0 Number  of  target  frequencies  in 

lining  (maximum  is  10).  If  NTFm 
is  set  equal  to  zero,  the  computer 
program  will  set  the  target  frequency, 
TFm,  to  the  current  calculated 
characteristic  frequency  for  a 
particular  noise  component.  The 
characteristics  frequency  is  that 
frequency  where  the  spectrum  level 
is  at  a maximum.  After  the 
default  target  frequency  is  set. 

NTFm  is  set  to  1 . 


NWLm 

0 

Number  of  walls  in  lining. 
(Maximum  is  10. > 

PCTAm(l) 

V 

100C 

Percent  treated  for  1st  target 
frequency. 

PCTAm(N) 

?( 

0 

Percent  treaed  for  Nth  target 
frequency,  N ^ NTFm. 

PLAm(  1 ) 

dB 

Peak  attenuation  for  1st  target 
frequency. 

PLAm(  N ) 

dB 

Peak  attenuation  for  Nth  target 
frequency,  N NTFm.  (Note:  Used 
only  if  IMAm  1 . ' 

R1  Wm(  1 ) 

m(ft) 

Radius  of  1 st  wall  of  lining. 

Rl  Wm(  N i 

m(ft) 

Radius  of  Nth  wall  of  lining. 
(Note:  Used  only  if  LGMm  = 1, 
N - NWLm.) 

TLm(  1 ) 

m(ft) 

Treatment  length  of  1st  wall  of 
lining,  N - NWLm. 

TLm)  N ) 

m(  ft) 

Treatment  length  of  Nth  wall. 

TFm(  1 ) 

Hz 

Charac- 
teristic 
peak  noise 
frequency 

First  target  frequency.  (Note:  If 
NTFm  0,  TFnU  1 ) w ill  he  reset  to  th 
current  calculated  characteristic 
frequency. ) 
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Variable 

Name 


Units 


Default  Description 


► 

f 

i 


TFmtNi 

Hz 

0 

Nth  target  frequency.  N = NTFm. 

Configuratio 

n Corrections: 

ICORm  f 0 

DOB(l) 

dB 

Table  of  AdB  corrections  as  a function 
of  directivity  angles  only.  (Note: 

DOB(N) 

dB 

ICORm  = 1 and  N = NPSCR). 

DPB(li 

dB 

Table  of  AdB  corrections  as  a function 
of  frequency  band  number,  and  the 

DPB(M) 

dB 

directivity  angles.  This  table  is 

input  as  a single  array  whose  indices 
correspond  to  a two-dimensional  array: 

DPB(  M)  = Xii.  j)  where  M = i + K (j-li,and 
(i,  j)  correspond  to  the  pass  band 
number  and  directivity  angle  index, 
respectively.  Note  that  K = 8 for  full 
octaves  or  24  for  1 3 octaves.  Applies 
only  if  ICORm  = 2. 


PSCR(l) 

deg 

Table  of  directivity  angles  corres- 
ponding to  either  DOB  array  or  DPB 
array  depending  on  IC'ORm  setting. 

PSCR(N) 

deg 

N NPSCR 

NPSCR 

Number  of  directivity  angles  on  which 
the  configuration  correction  table  is 
based.  (2  - NPSCR  - 17*. 

4.3.2  PRIMARY  JET  NOISE  AND  PRIMARY  PLUS  SECONDARY  JET  NOISE 

This  section  describes  the  subset  of  the  SNOISIN  parameters  used  in  estimating  either 
primary  jet  noise  or  combined  primary  and  secondary  jet  noise  These  inputs  are  needed 
in  addition  to  the  appropriate  $GDATA  parameters  described  previously  (Note:  For 
AdB  corrections,  see  sec  4 2 1 on  lining  attenuation  and  configuration  corrections.) 
Applicable  shielding  data  for  this  subset  are  described  in  section  4.4 
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Variable 

Name 


1 ' I ! I t .- 


I letaillt 


I )e-.i  ri|.i  mi 


’ INSHLD  n 


ITYPE 


NTYPE  1 


Primary  Jet  Parameters: 

API  m2(ft2) 

ANGJTl  deg  0 

DIAMTl  m(ft)  0 


N.JETl 


I mill  al.  ir  (■  >r  del  ml  mg  wing  .shielding 
Tab  ulations  are  desired  and  inputs 
are  required  . mrnediatel v after  this 

SNOISIN  data  set 

INSHLD  (t  for  nu  input  or  shielding 
calculation 

INSHLD  / 0 for  input  according  to 

$SHLDAT  data  set  format 
and  shielding  calculations 
are  desired 

RESTRICTION:  If  default  values  are 
desired  in  shielding  calcualtions. 

INSHLD  - 1 must  be  input  on  the  first 
$N()ISIN  data  set  in  the  first  case 
of  a job. 

Indicator  for  primary  or  combined 
primary  and  secondary  jet  noise. 

ITYPE  --  1 for  primary  jet 

ITYPE  = 2 for  primary  and  secondary  jet 

This  variable  must  be  specified  in  the 
first  case  for  each  noise  component 
for  each  different  noise  source 
configuration. 

Number  of  noise  components  in  a 
configuration.  (Note:  NTYPE  must 
be  specified  only  in  the  first 
SNOISIN  data  set  of  each  configuration 
in  the  first  case  of  a run  if 
different  than  1 . ) 


Cross-sectional  area  of  the  nozzle 
exit.. 

Engine  inclination  angle. 

Diameter  of  nozzle. 

If  zero  or  negative,  the  diameter 
will  be  calculated  based  on  nozzle 
area  (API  >. 

Code  for  type  of  input  data. 

NJET1  = 1 user  defines  API , PRl , TT1 
N.JETl  2 user  defines  WP1.  PRL  TT1 
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Variable 

Name 


Units 


Default  Description 


- 

NJETl  = 3 user  defines  API,  WP1,  VP1 
plus  AS2,  VS2,  WS2  if 
secondary  jet  noise  is 
to  be  considered. 

MCODE1 

1 

Code  for  Strouhal  curve. 

MCODE1  = 1 for  flight  spectrum  curve 
MCODE1  = 2 for  ground  spectrum  curve. 

PR1 

Nozzle  pressure  ratio;  i.e.,  total 
pressure  divided  by  freestream 
static  pressure. 

TT1 

K(°R) 

Jet  total  temperature. 

VP1 

m/s(fps) 

Velocity  of  jet  exhaust  relative  to 
nozzle. 

WP1 

kg/s(lbm/s) 

Primary  mass  flow. 

Secondary  Jet  Parameters: 

The  following  three  parameters  are  needed  in  addition  to  the  above  combined  primary 
and  secondary  jet  noise,  (ITYPE  = 2). 


Variable 

Name 

Units  Default 

Description 

AS2 

m2(ft2) 

Secondary  jet  nozzle  area. 

VS2 

m/s(fp8) 

Secondary  jet  velocity  relative  to 
nozzle. 

WS2 

kg/sdbm/s) 

Secondary  mass  flow. 

4.3.3  CORE  AND  TURBINE  NOISE 


This  section  describes  the  subset  of  the  $NOISIN  parameters  used  to  estimate  core  and 
turbine  noise.  These  inputs  are  specified  in  addition  to  the  $GDATA  data  set  described 
previously.  (Note:  For  AdB  corrections,  see  sec.  4.3.1  on  lining  attenuation  and 
configuration  corrections.) 


Variable 

Name 

Units 

Default 

Description 

DELT3 

deg 

0 

Engine  attitude  angle. 

INSHLD 

0 

Wing  shielding  indicator.  Description 
same  as  ITYPE  = 1,  2,  3.  or  4. 

ISW3 

0 

Specifies  noise  type  to  be  predicted. 
ISW3  = 0 for  core  and  turbine  noise 

Variable 

Name 


Units 


Default  Description 


[ 

* 

i 

r 


ISW3  = 2 for  core  noise  only 
ISW3  = 3 for  turbine  noise  only 
RESTRICTION:  When  wing  shielding 
is  desired  (INSHLD  ^ 0),  the  value  of 
ISW3  will  impact  the  empirical  terms 
chosen  for  shielding  calculations. 

ISW3  = 2 causes  program  to  use  core 
noise  empericisms 

ISW3  ^ 2 causes  program  to  use  turbine 
noise  empericisms 


ITYPE 

ITYPE  = 3 for  core  and  turbine 
noise  prediction 

This  variable  must  be  specified  in 
the  first  case  for  each  noise 
component  for  each  configuration. 

NENG 

1 

Number  of  engines.  If  other’ than 
1 , this  must  be  specified  for  the  first 
noise  component  of  each  type  of 
propulsion  system. 

NTYPE 

1 

Number  of  noise  types  in  a configure 
tion.  (Note:  NTYPE  must  be  speci- 
fied only  in  the  first  $NOISIN  data 
set  of  each  configuration  for  the 
first  case  of  a run.) 

Core  Noise  Parameters: 

CMF3 

kg/s(lbm/s) 

Combustor  corrected  mass  flow. 
Corrected  to  sea  level,  static 
conditions  (1  ATM,  15°  C). 

EK3 

Specific  engine  correction.  (See 
table  10  in  NASA  CR1 14649.) 

JB3 

Indicator  for  type  of  burner. 

JB3  = 1 for  annular  type  burner 
JB3  = 2 for  can  type  burner 

PP3 

Turbine  total  pressure  ratio;  i.e., 
turbine  inlet  total  pressure  divided 
by  turbine  exit  total  pressure. 

TT3 

K(°R) 

Combustor  exit  total  temperature. 

1 
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Variable  Units  Default  Description 
Name 

Turbine  Parameters: 


BN3 

Number  of  blades  for  turbine  last 
stage. 

CLS3 

m/s(fps) 

Speed  of  sound  at  last  turbine  stage 
If  CLS3  is  not  set  it  will  be 
estimated  internally  by  program. 

CS3 

Stator/rotor  spacing.  (See  fig. 
55  in  NASA-CR1 14649.) 

DT3 

m(ft) 

Tip  diameter  for  turbine  last  stage 
rotor.  Required  only  if  VTR3  is 
unknown. 

IC3 

Indicator  for  nozzle  configuration 
type. 

IC3  = 0 for  dual  flow  nozzles  of 
turbofans  or  turbojets 
IC3  ? 0 for  engines  with  retracted 
primary  flow  nozzle 
(e.g.,  JT8D > 

PMF3 

kg/8(lbm(s) 

Primary  mass  flow. 

SS3 

RPM 

Shaft  speed. 

TU3 

K(°R) 

Turbine  outlet  total  temperature. 
Required  only  if  CLS3  is  unknown. 

VTR3 

m/s(fps) 

Relative  tip  speed  of  turbine  last 
stage  rotor. 

4.3.4  COMPRESSOR,  FAN  INLET,  AND  FAN  EXIT  NOISE 

This  section  describes  the  subset  of  the  $NOISIN  parameters  used  to  estimate  the 
compressor,  fan  inlet,  and  fan  exit  noise.  These  inputs  are  needed  in  addition  to  the 
appropriate  $GDATA  parameters.  (Note:  For  AdB  corrections,  see  sec.  4.3.1  on  lining 
attenuation  and  configuration  corrections.) 


Variable 

Name 

Units 

Default 

Description 

DELT45 

deg 

0 

Engine  attitude  angle. 

FPR45(I) 

Fan  or  compressor  pressure  ratio, 
(1  s I s£  NSTG45). 

INSHLD 

0 

Wing  shielding  indictator.  Descript- 
ion is  same  for  1TYPE  = 1 , 2,  3.  or  4 
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Variable 

Name 


Units 


Default  Description 


I 


1 


ITYPE 


NB45(I) 

NENG  1 


NTYPE  1 


NSTG45 

RN145  RPM 

RSS45(I)  % 

RTS45  0 


Fan  Inlet  Parameters: 
CFPR4 

DIAM4(I)  m(ft) 


Fan  Exit  Parameters: 


AREA5(I)  m2(ft2) 
BPR5 


ITYPE  = 4 for  compressor  or  fan 
inlet  noise 

ITYPE  = 5 for  fan  exit  noise 

This  variable  must  be  specified  in 
the  first  case  for  each  noise 
component  for  each  configuration. 

Number  of  compressor  or  fan  blades 
for  each  stage,  (1  I s=  NSTG45). 

Number  of  engines.  If  other  than  1,  this 
must  be  specified  for  the  first 
noise  component  of  each  type  of 
propulsion  system. 

Number  of  noise  components  in  the 
configuration.  (Note:  NTYPE  must 
be  specified  only  in  the  first 
$NOISIN  data  set  of  each  configuration 
for  the  first  case  of  a run.) 

Number  of  fan  stages,  (1  « NSTG45  « 3). 
Rotor  rotational  speed. 

Minimum  rotor/stator  spacing, 

(1  « I NSTG45). 

Relative  tip  Mach  number  of  the 
first  stage  without  inlet  guide 
vanes  (IGV).  If  less  than  or 
equal  to  0,  IGV’s  will  be  assumed 
for  the  first  stage.  If  RTS45  is 
less  than  one,  but  greater  than 
zero,  there  is  no  buzzsaw  component. 

(In  addition  to  the  inputs  above  for 
inlet  fan  noise;  i.e.,  ITYPE  = 4.) 

Fan  pressure  ratio  when  the  relative 
tip  Mach  number  equals  1.025. 

Compressor  or  inlet  fan  diameter, 

(1  s I « NSTG45). 

(In  addition  to  the  first  set  of 
inputs  for  fan  exit  noise;  i.e., 

ITYPE  = 5.) 

Fan  discharge  area,  (Isis  NSTG45). 

Engine  bypass  ratio,  m2/m|  where 
mj  . . . primary  mass  flow 
m2  . . . secondary  mass  flow 
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Variable  Units  Default  Description 
Name 

Indicator  for  duct  type 
= 0 for  short  fan  ducts 
= 1 for  long  fan  ducts  with  retracted 
primary  nozzle;  i.e.,  the  JT8D 
engine 

= 2 for  long  fan  ducts  with  approx- 
imate coplanar  primary /secondary 
nozzle  exits 

= 3 for  approximate  3/4  length  fan 
ducts 


4.3.5  AUGMENTOR  WING  NOISE 

This  section  describes  the  subset  of  the  $NOISIN  parameters  used  to  estimate  the 
augmentor  wing  noise.  These  inputs  are  needed  in  addition  to  the  appropriate  $GDATA 
parameters  described  previously.  (Note:  For  AdB  corrections,  see  sec.  4.3.1  on  lining 
attenuation  and  configuration  corrections.) 


Variable 

Name 

Units 

Default 

Description 

AD6 

m2(ft2) 

Nozzle  discharge  area. 

DE6 

m(ft) 

Effective  diameter  (hydraulic  diameter) 
DE6  = 4 * Al)6  perimeter 
= 2 H <1  + H L) 

H is  slot  height;  L is  slot  length. 

DELT6 

deg 

0 

Flap  angle  relative  to  the  horizon. 

GAMA6 

1.4 

Ratio  of  specific  heats  for  exhaust 
flow. 

ITYPE 

ITYPE  ~ 6 for  augmentor  wing  noise. 
This  variable  must  be  specified  in 
the  first  case  for  each  noise 
component  for  each  configuration. 

NENG 

1 

Number  of  engines.  If  other  than  1,  this 
must  be  specified  for  the  first  noise 
component  of  each  type  of  propulsion 
system. 

NTYPE 

1 

Number  of  noise  components  in  a 
configuration.  (Note.  NTYPE  must  be 
specified  only  in  the  first  $NOISIN 
data  set  of  each  configuration  for 
the  first  case  of  a run.) 

4 

i 


39 


/ 


i 

i 

I 


i 


< 


i 

i 


< 


/ 

I 


i 


Variable 

Name 

Units 

Default 

Description 

TT6 

K(°R) 

Total  temperature  at  the  nozzle  exit. 

XNPR6 

Nozzle  pressure  ratio-the  total 
pressure  at  the  nozzle  exit  divided  by 
the  freestream  static  pressure. 

4.3.6  BLOWN  FLAP  NOISE 

This  section  describes  the  subset  of  the  SNOISIN  parameters  used  to  estimate  the  blown 

flap  noise. 

These  inputs 

are  supplied 

in  addition  to  the  $GDATA  parameters.  (Note:  For 

AdB  corrections,  see  sec. 

4.3.1  on  lining  attenuation  and  configuration  corrections.) 

Variable 

Name 

Units 

Default 

Description 

AN7 

ml(ft2) 

Nozzle  discharge  area. 

DELT7 

deg 

0 

Engine  attitude  angle. 

DL7 

Dimensionless  distance  between 
nozzle  exit  and  target  point  on  the 
flap(s)  when  the  nominal  flap  angle 
is  45°;  i.e.,  L/D  in  NASA-CR1 14649 
(sec.  5.2.2.51. 

DN7 

m(ft) 

Nozzle  exit  diameter  or  hydraulic 
diameter. 

FANG7 

deg 

0 

Nominal  flap  angle. 

HD7 

Dimensionless  distance  between  nozzle 

centerline  and  mean  wing  chord;  i.e., 
H/D  in  NASA-CRl  14649  (sec.  5.2.2.5). 

1TYPE 

ITYPE  = 7 for  blown  flap  noise.  This 
variable  must  be  specified  in  the  first 
case  for  each  noise  type  for  each 
configuration. 

NENG 

1 

Number  of  identical  noise  sources. 

If  other  than  1,  this  must  be  specified 
for  the  first  case  for  each 
configuration. 

NTYPE 

1 

Number  of  noise  types  in  a configura- 

tion.  (Note:  NTYPE  must  be  specified 
only  in  the  first  $NOISIN  data  set  of 
each  configuration  for  the  first  case 
of  a run.) 
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Variable 

Name 

Units 

Default 

Description 

PR7 

- 

Nozzle  pressure  ratio;  i.e.,  total 
pressure  divided  by  freestream  static 

pressure. 

TT7 

K(°R) 

Total  temperature  of  exhaust  at 
nozzle  exit. 

4.3.7  LIFT  FAN  NOISE 

This  section  describes  the  subset  of  the  $NOISIN  parameters  used  to  estimate  lift  fan 
noise.  These  inputs  are  supplied  in  addition  to  the  appropriate  SGDATA  parameters 
described  previously.  (Note:  For  AdB  corrections,  see  sec.  4.3.1  on  lining  attenuation  and 
configuration  corrections.) 


Variable 

Name 

Units 

Default 

Description 

AREA8 

mW) 

Fan  discharge  area.  IF  = 0.  No  aft 
fan  noise  is  calculated. 

CRFPR8 

Fan  pressure  ratio  for  the  relative 
tip  Mach  number  of  1.025. 

DELTA8 

deg 

0 

Engine  attitude  angle. 

DIAM8 

m(ft) 

Fan  inlet  diameter.  IF  = 0.  No  inlet 
fan  noise  calculated. 

FPR8 

Fan  pressure  ratio;  i.e.,  total 
pressure  aft  of  a fan  stage 
divided  by  total  pressure  just 
forward  of  the  fan  stage. 

ITYPE 

ITYPE  = 8 for  lift  fan  noise.  This 
variable  must  be  specified  in  the 
first  case  for  each  noise  type  for 
each  configuration. 

NB8 

Number  of  fan  blades. 

NENG 

1 

Number  of  lift  fans  being 
considered.  If  other  than  1 , this  must 
be  specified  for  the  first  ease  for 
each  configuration. 

NTYPE 

1 

Number  of  noise  types  in  a configura- 
tion. (Note:  Must  be  specified  only 
in  the  first  SNOISIN  data  set  of  each 
configuration  for  the  first  case  of  a 
run.) 

RN18 

RPM 

Rotor  rotational  speed. 
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Variable 

Name 

Units 

Default 

Description 

RSS8 

% 

Minimum  rotor/stator  spacing. 

RTS8 

0 

Relative  tip  Mach  number  of  the  fan 
without  inlet  guide  vanes.  If  RTS8 
is  less  than  or  equal  to  zero  inlet 
guide  vanes  are  assumed.  If  less 
than  one  but  greater  than  zero, 
there  is  no  buzzsaw  noise  component 

4.3.8  EJECTOR-SUPPRESSOR  NOISE 


This  section  describes  the  subset  of  the  $NOISIN  data  set  used  to  estimate 
ejector-suppressor  noise.  These  inputs  are  supplied  in  addition  to  the  $GDATA 
parameters  described  previously,  (Note:  For  AdB  corrections,  see  sec.  4.3.1  on  lining 


attenuation  and  configuration.) 

Variable 

Name 

Units 

Default 

Description 

AR9 

Area  ratio  of  suppressor  nozzle;  i.e., 
primary  plus  secondary  flow  area 
divided  by  primary  flow  area. 

AREA9 

m2(ft2) 

Discharge  area  of  suppressor  nozzle. 

CV9 

Velocity  coefficient  for  nozzle. 

DELT9 

deg 

0 

Angle  between  thrust  vector  and 
horizon. 

EMACH9 

Exhaust  Mach  number  for  ejector 
(only  needed  if  IEJ9  ^ 0). 

EXNM9 

Exhaust  Mach  number  for  nozzle. 

ITYPE 

ITYPE  = 9 for  ejector-suppressor  noise. 
This  variable  must  be  specified  in 
the  first  case  for  each  noise  type 
for  each  configuration. 

IEJ9 

0 

Switch  for  ejector  and/or  suppressor. 
IEJ9  = 0 bare  suppressor 
IEJ9  ^ 0 ejector/suppressor 

NENG 

1 

Number  of  engines.  If  other  than  1,  this 
must  be  specified  for  the  first  case 
for  each  configuration. 

NTYPE 

1 

Number  of  noise  types  in  a configura- 
tion. (Note:  NTYPE  must  be  specified 
only  in  the  first  $NOISIN  data  set  of 
each  configuration  for  the  first  case 
of  a run.) 
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Variable 

Name 

Units  Default 

Description 

NUM9 

- 

Number  of  discharge  elements  of 
suppressor  nozzle. 

PA9 

m2(ft2) 

Discharge  area  of  ejector  (required 
only  if  IEJ9  / 0). 

PCV9 

Velocity  coefficient  for  ejector 
(required  only  if  IEJ9  ^ 0). 

PS9 

ATM(psia) 

Static  pressure  in  exhaust  at 
nozzle  (required  only  ll  1EJ9  / 0). 

SMACH9 

Mach  number  of  induced  secondary  air 

ST9 

K(°R) 

Static  temperature  at  nozzle  exit. 

4.3.9  PROPELLER  NOISE 


This  section  describes  the  subset  of  the  $NOISIN  parameters  used  to  estimate  propeller 
noise  using  the  emprical  procedure  defined  in  NASA-CRl  14649.  The  next  section 
describes  inputs  for  the  theoretical  rotor  procedure  which  may  be  used  in  lieu  of  this 
module.  These  inputs  are  needed  in  addition  to  the  appropriate  $GDATA  parameters 
described  previously.  (Note:  For  AdB  corrections,  see  sec.  4.3.1  on  lining  attenuation  and 
configuration  corrections. ' 


Variable  Units  Default  Description 

Name 


ASUBIO  m2(ft2) 

BIO 

DIO  m(ft) 

DELT10  deg  0 

DSUB10  m(ft) 

ITYPE 


NENG  1 

NTYPE  1 


Total  blade  area  for  one  side  of 
propeller. 

Number  of  propeller  blades. 
Propeller  diameter. 

Angle  between  gross  thrust  vector 
and  horizon. 


Characteristic  dimension  for  the 
blade  geometry  at  0.7  span;  i.e., 
the  axial  projected  chord. 

ITYPE  --  10  for  propeller  noise.  This 
variable  must  be  specified  in  the 
first  case  for  each  noise  type  for 
each  configuration. 


Number  of  engines.  If  other  than  one,  this 
must  be  specified  for  the  first  case 
for  each  configuration. 


Number  of  noise  types  in  a given 
configuration.  (Note:  NTYPE  must 
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Variable 

Name 


Units 


Default  Description 


- 

be  specified  only  in  the  first 
$NOISIN  data  set  of  each  configuration 
for  the  first  case  of  a run). 

RPM10 

RPM 

Propeller  rotational  speed. 

T10 

N(lbf) 

Propeller  thrust. 

W10 

kW(HP) 

Propeller  shaft  power. 

4.3.10  HELICOPTER,  PROPELLER,  AND  TILT  ROTOR  NOISE 

This  section  describes  the  subset  of  the  $NOISIN  parameters  used  to  estimate 
helicopter,  propeller,  and  tilt  rotor  noise  based  on  the  theoretical  procedures  defined  in 
section  5. 2. 5. 2 of  NASA  CR1 14649.  These  inputs  are  needed  in  addition  to  the 
appropriate  $GDATA  parameters  described  previously.  (Note:  For  AdB  corrections,  see 
sec.  4.3.1  on  lining  attenuation  and  configuration  corrections.) 


Variable 

Name 

Units 

Default 

Description 

AB11 

Bll 

m2(ft2) 

Total  blade  area  of  one  side  of  rotor. 

Number  of  blades  per  rotor, 

(2.  s Bll  =£  6.). 

CEE11 

24.4 

Constant  (c)  in  loading  law  (eq.  51  of 
NASA-CR114649).  This  variable  and 
variables  XMMll  and  XLMC11  must 
be  specified  when  LLF11  = 6. 

DELTll 

deg 

Angle  between  gross  thrust  vector 
and  horizon.  (Applies  only  to 
the  main  rotor.) 

DE11 

m(ft) 

Characteristic  dimension  for  the  blade 
geometry  at  0.7  span;  i.e.,  the  mean 
axial  projected  chord. 

DT11 

ITYPE 

m(ft) 

Tip  diameter. 

ITYPE  =11  for  helicopter  noise.  This 
variable  must  be  specified  in  the 
first  case  for  each  noise  type  for 
each  configuration. 

IRR11 

0 

Indicator  for  specifying  if  the  rotor 
being  considered  is  the  main  rotor  or 
tail  rotor. 

IRR11  = 0 for  main  rotor 
IRR1 1 ^ 0 for  tail  rotor 
(Note:  If  the  tail  rotor  is  being 

rj*.- 
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Variable 

Name 


Units  Default  Description 


LLF11 


considered,  the  thrust  axis  is  assumed 
horizontal  and  perpendicular  to  the 
helicopter's  flightpath.) 

Loading  law  indicator.  (See 
NASA  CR2 14649  for  equations.) 

LFll  = 1 applying  to  hovering 
helicopter  (eq.  50A) 

LFll  = 2 applying  to  helicopters 
and  tilt  rotors  (eq.  50B) 
LFll  = 3 applying  to  low  speed 
propellers  (eq.  50C) 

LLF11  = 4 applying  to  low  speed 
propellers  (eq.  501)  > 

LLF11  = 5 applying  to  medium  speed 
propellers  (eq  5uE> 

LFll  = 6 user  inputs  loading  law 
parameters  (eq.  51 ) 


NENG 

1 

Number  of  engines;  i.e..  rotors.  If 
other  than  1 , this  must  be  specified 
for  the  first  noise  component  of 
each  type  of  propulsion  system. 

NTYPE 

1 

Number  of  noise  types  in  the  configura 
tion.  (Note:  NTYPE  must  be  specified 
only  in  the  first  SNOISIN  data  set  of 
each  configuration  for  the  first  case 
of  a run.) 

Qll 

N-m(ft-lbf) 

Shaft  torque 

RN11 

0.8 

Dimensionless  centroid  for  equivalent 
point  load  on  a rotor  blade. 

RPM11 

RPM 

Rotor  rotational  speed 

Sill 

5.0 

Lift  curve  slope  for  a single  blade 
'applies  if  I.l.Fl  1 2i 

Til 

N(lbf) 

Thrust  Her  rotor 

XMM11 

2. 

• '(instants  tin  & A,.'  in  loading  law 
(eq  51  of  NASA  CR  1 14649). 

XLMC11 

30.  . 

Applies  when  LLFl  1 6. 

4.3.11  MEASURED  DATA  INPUT 


This  section  describes  the  subset  of  the  $NOISlN  parameters  required  for  inclusion  of 
measured  data.  The  SPL  variable  described  in  this  section  is  an  array  of  sound  pressure 
levels  in  dB  re  20  /xN/'m2  as  a function  of  frequency  (preferred  1 1 octave  bands  or  1 3 


45 


octave  bands),  a prescribed  engine  performance  parameter,  directivity  angle,  and 

elevation  angle.  In  order  to  minimize  core  size,  this  SPL  array  is  assigned  to  local  . 

storage  in  the  measured  data  overlay  and  the  SPL  input  is  read  in  an  eight  field,  10 
column  per  field  decimal  format  rather  than  the  NAMELIST  $NOISIN  format.  The 

measured  data  discussed  in  this  section  is  limited  to  one  powerplant  and  can  be  ^ 

specified  only  once  during  a computer  run.  The  only  variables  which  can  change  in  i 

succeeding  cases  are  the  $GDATA  parameters,  DELT12,  and  the  configuration  \ 

corrections  discussed  in  section  4.3.1.  \ 

j 

Refer  to  section  3.3.2  of  NASA-CR114650  for  sample  input  data  which  includes  j 

measured  data.  9 

Variable  Units  Default  Description  j 

Name  | 


BETA12(1) 

deg 

Independent  variable  array  of 
elevation  angles  used  to  correlate 

BETA12(N) 

deg 

the  measured  SPL’s.  Axial  symmetric 
sound  sources  will  have  no  BETA12 
entries  indicated  by  NBTA12  = 0. 

N = NBTA12. 

DELT12 

deg 

0 

Engine  inclination  angle. 

EP12(1) 

Independent  variable  array  of  engine 
perforamnce  parameters  needed  to 

EP12(N) 

correlate  the  measured  SPL’s.  N = NEP12 

ITYPE 

ITYPE  = 12  for  measured  data.  This 
variable  must  be  specified  in  the 
first  case  for  each  noise  type  for 
each  configuration. 

NBTA12 

0 

Number  of  enties  in  the  BETA12  data 
array,  (NBTA12  = 0 or  2 « NBTA12  s 5). 

NEP12 

Number  of  entries  in  the  EP12  data 
array,  (2  =£  NEP12  =s  5). 

NPSI12 

Number  of  entries  in  the  PSI  12  data 
array,  (2  =s  NPSI12  s 17). 

NTYPE 

1 

Number  of  noise  types  in  a configura- 
tion. (Note:  NTYPE  must  be  specified 
only  in  the  first  $NOISIN  data  set  of 
each  configuration  for  the  first  case.) 

PSI  12(1) 

deg 

Independent  variable  array  of 
directivity  angles  used  to  correlate 

PSI12(N) 

deg 

the  measured  SPL’s.  N = NPSI12. 

(Note:  SPL  array  is  input  with  8F10.0 
format  directly  following  the 
$NOISIN  data  set.) 
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Variable 

Name 


Units 


Default  Description 


\ 


E 


1 


f 


Dependent  data  array  of  SPL’s  in  dB 
for  noise  versus  (f.  PSI12,  EP12, 

BETA12).  Where  f is  the  frequency 
(eight  preferred  1/1  octave  bands  or 
twenty-four  1/3-octave  bands)  defined 
in  table  4 of  NASA-CR1 14649.  Inputs 
are  for  free-field  noise  conditions 
at  R = 1 m. 

Note:  SPL(n)  = F(I,  J,  K,  L)  where 
n — I -f-  kj  (*J-1  -f-  k2  (K-l  + kt3 ) (L-l)t 
kj  = 24  for  1/3-octave  band  analysis 
= 8 for  full  octave  band  analysis 
ko  = NPSI12 
k3  =NEP12 

(I,  J,  K,  L)  are  indices  corresponding  to 
(f,  PSI12,  EP12,  BETA12)  where  frequency 
is  varied  first,  then  the  directivity 
angle,  then  the  engine  performance  para- 
meter, and  finally  the  elevation  angle. 

4.3.12  JET/EDGE  INTERACTION  NOISE* 

This  section  describes  the  subset  of  the  $NOISIN  data  set  used  to  estimate  jet  edge 
interaction  noise.  These  inputs  are  needed  in  addition  to  the  appropriate  SGDATA 
parameters  described  previously.  (For  AdB  correction,  refer  to  sec.  4.3.1  lining 
attenuation  and  configuration  corrections.) 


Variable 

Name 

Units 

Default 

Description 

ITYPE 

ITYPE  =13  for  jet  edge  interaction 
noise  prediction.  This  variable 
must  be  specified  in  the  first 
case  for  each  noise  component  for 
each  configuration. 

NENG 

1 

Number  of  engines.  If  other  than  1,  this 
must  be  specified  for  the  first 
noise  component  of  each  type  of 
propulsion  system. 

NTYPE 

1 

Number  of  noise  components  in  the 

configuration.  (Note:  NTYPE  must 
be  specified  only  in  the  first  NOISIN' 
data  set  of  each  configuration  for 
the  first  case  of  a run.) 


SPUl) 

dB 

SPUN) 

dB 

1 


47 


1 I 


i>  ' 


Variable 

Name 

EDGVAR(l) 


Units 


Default  Description 


Nondimensional  distance  between 
nozzle  exit  to  the  wing/flap  systems 
trailing  edge  measured  along  the 
wing/flap  surface  (L/D). 

Jet  Mach  number. 

Hydraulic  diameter  of  jet  nozzle  (D) 
(equals  four  times  discharge  area 
divided  by  nozzle  perimeter). 

Jet  static  temperature  ratio 
(Ts/Ts0  whsre  Tso  is  288.2  K or 
518.7  °R). 

Not  used. 

Jet  spreading  angle. 

Nondimensional  discharge  area  (A/D2). 

Nondimensional  nozzle  lip  height 
above  wing  surface  (H/D). 

Nondimensional  jet  core  length  (Lc/D). 
The  default  is  for  round  nozzle. 

Nominal  flap  angle  or  angle  of 
attack  when  the  flaps  are  retracted. 


4.4  NOISE  COMPONENT  SHIELDING  PARAMETERS* 

($SHLDAT  INPUT  DATA  SET) 

If  noise  shielding  is  specified  for  the  particular  engine/wing  configuration  (INSEOW  in 
$GDATA)  and  for  the  particular  noise  components  (INSHLD  = 1 in  $NOISIN)  then  the 
data  set  $SHLDAT  must  follow  the  INOISIN  data  set.  The  $SHLDAT  data  set  contains 
the  shielding  inputs  for  the  noise  component  specified  in  the  $NOISIN  data  set.  The 
noise  components  which  can  include  shielding  effects  are: 


1 

EDGVAR(2) 

I 

EDGVAR(3) 

m(ft) 

EDGVAR(4) 

l ' 

r 

\ 

EDGVAR(5) 

EDGVAR(6) 

deg 

5.7 

EDGVAR(7) 

0.7854 

r. 

f 

EDGVAR(8) 

0. 

[ 

l 

EDGVAR(9) 

5. 

1 

l . 

\ 1 

EDGVAR(IO) 

deg 

0. 

Noise  Component 

Primary  jet 

Primary  and  secondary  jtt 
Core  and  turbine 
Inlet  fan  or  compressor 
Exit  fan 


ITYPE 

1 

2 

3 

4 

5 
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Variable 

Name 


Units 


Component  Default  Description 
ITYPE 


l< 


t 

\ 

r 

i 


1 

i ‘ 


i i 
1 1 


ASF  deg  3,5 


BETA(Ii  3,5 


CPSIO(I)  3,5 


DIAMT2  m(ft)  2 


DSL1 


1,2 


13.  Empirical  constant  y 

in  formula  for  computing 
the  angular  shift  of 
the  apparent  sound 
source  due  to  exhaust 
flow.  (See  Fig.  A-10.) 

Fig.  A-10  One-dimensional  array 
used  to  define  the 
dependent  variable  fi 
in  the  emprical  curve 
fi  versus  cotdji,,)  of 
Figure  A-10.  Values  of 
BETA  represent  the 
radial  offset  of  a 
hypothetical  line 
parallel  to  the  wing’s 
equivalent  half-plane 
used  to  locate  the 
apparent  sound  source, 

(1  « I s INASRO). 

Fig.  A-10  One-dimensional  array 
used  to  define  the  in- 
dependent variable 
cotliji,,)  in  the  emprical 
curves  f3  versus  cot(i/<0)  and 
a versus  cot(ifi()).  Values 
of  CPSIO  represent  the 
cotangent  of  the  direct- 
ivity angle  rela- 
tive to  the  engine 
inlet  centerline,  graz- 
ing the  wing  edge  for  the 
shortest  sound  propa- 
gation path  from  the  nozzle 
exit  to  the  edge  to  the 
observer,  (2  «s  I « INASRO). 

Diameter  of  nozzle 
(secondary  jet). 

If  noncircular,  use 
hydraulic  diameter. 

Nondimensional  shield  length 
measured  parallel  to  the  exhaust 
axis  from  the  nozzle  exit  plane. 
(Lj  D i ) primary  jet. 


i 


] 


$ 


a 
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FLYOVER  ILLUSTRATION 
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Figure  AW . - Empirical  Adjustments  to  Jet  Refraction  Theory  for 
Apparent  Source  Position 


Variable 

Name 


Units 


Component  Default  Description 

ITYPE 


1 

I 

I 

I 


i 

t- 


DSL2 

.2 

Nondimensional  shield 
length  measured  parallel 
to  the  exhaust  axis 
from  the  nozzle  exit 
plane. 

(L2/D2)  secondary  jet. 

EMJ 

3,5 

Exhaust  flow  Mach  number 

FASS(I) 

3,5 

Fig.  A-10  One-dimensional  array 
used  to  define  the 

dependent  variable 

in  the  empirical  curve 

a versus  cotliji,,)  of 

figure  A-10.  Values  of 

FASS  represent  the 

correction  fautor  (a) 

in  the  formula  for 

computing  the  angular 

shift  of  the  apparent 

sound  source  due  to 

exhaust  flow,  (2  I INASRO) 

INASRO  3,5  0 Number  of  entries  in 

the  input  arrays  BETA, 

CPSIO,  and  FASS  used 

to  define  the  empirical 

curves  (3  versus  cot(tJ/0) 

and  a versus  cot(<i/n)  similar 

to  that  shown  in 

figure  A-10.  If  INASRO  is 

not  input,  then  the 

default  value  0 signals 

the  computer  to 

initialize  the  values 

for  BETA,  CPSIO.  and 

FASS  to  correspond  to 

figure  A-10,  (2  t INASRO  *=  24 

INUSP  3,4,5  0 Indicator  to  denote 

which  empirical  direct- 
ivity curve  to  be  used 
for  shielding  calculations. 

= 0 for  predicted 
USPL  versus  i h 
- 1 for  built-in  curves 
of  figure  A- 11. 

= 2 for  user  input  curves 
USPL  versus  PSI 
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Variable 

Name 

Units 

Component 

ITYPE 

Default 

IWED(l) 

3,4,5 

0 

IWED(2) 

0 

IWED(3) 

0 

NUSPL 

3,4,5 

See  note 
below 

PSI(I) 

deg 

3,4,5 

See  note 
below 

Description  j 

| 

Array  of  values  equal  3 

to  0 or  1 where  1 

indicates  that  the  wing  edge  is  i 

to  be  considered  in 
shielding  calculations. 

IWED(  1 ) denotes  trailing  edge  i 

IWED12)  denotes  leading  edge  ‘ 

IWED(3)  denotes  tip  edge 

The  number  of  entires  in 
the  empirical  directivity 
curve  USPL  versus  PSI. 

3 « NUSPL  =s  19. 

An  array  of  directivity 
angles  in  the  empirical 
directivity  curve 
USPL  versus  PSI, 

(UU  NUSPL1. 


1,  i 

TSTSO 

3,5 

1.  Exhaust  flow  static 

; 

temperatue  ratio 

< TS/TS(>  i where  Tsn  is  the 

f , 

i 1 

ambient  static  tempera- 

1 

. ) 

ture  in  absolute  units 

: i 

(°R  or  Kl. 

4 


USPL(I)  dB  3,4,5 


See  note  An  array  of  sound 

below  pressure  levels  for  the 

empirical  unshielded 
directivity  curve 
USPL  versus  PSI, 


(1  « I s,  NUSPL). 

Note:  The  default  values  of  the  curve  USPL  versus  PSI  for  the  various  components  are 
based  on  static  tests.  (See  fig.  A-ll.) 


1 


1 
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5.0  OUTPUT  DATA 


Output  data  are  recorded  on  the  OUTPUT  file  after  the  execution  of  each  case,  that  is, 
for  each  new  set  of  aircraft  conditions.  The  user  has  the  option  of  requesting  up  to  seven 
different  reports  as  described  in  section  4.1  under  General  Data  Parameters.  The 
following  sequence  of  output  is  printed  per  case  when  all  reports  are  requested. 

1.  Total  noise  spectra  (all  components)  observed  at  each  sideline  distance. 

2.  Assumptions  for  noise  estimation,  a summary  of  options. 

3.  First  noise  component  spectra  observed  at  each  sideline  distance  for  the  first 
engine  configuration.  Second  noise  component  spectra  observed  at  each  sideline 
distance  for  the  first  configuration. 

Nth  noise  component  spectra  observed  at  each  sideline  distance  for  the  first  engine 
configuration. 

Nth  noise  component  spectra  observed  at  each  sideline  distance  for  the  next  engine 
configuration. 

4.  Flightpath/observer  geometry 
engine/wing  geometry  (first  configuration). 

engine/wing  geometry  (next  configuration,  etc.) 

5.  Noise  extrapolation  corrections. 

6.  Total  noise  spectra  at  index  conditions  emitted  toward  the  observer  position  at 
each  sideline  distance. 

7a.  First  noise  component  spectra  at  index  conditions  emitted  toward  the  observer 
positioned  at  each  sideline  distance  for  the  first  engine  configuration. 

7b.  If  applicable,  estimated  configuration  corrections  which  can  include  measured 
corrections,  lining  corrections  and  wing  shielding  corrections. 

7c.  Same  except  for  the  next  noise  component  of  that  engine  configuration. 

Repeat  of  7a  through  7c  except  for  the  next  engine  configuration. 

This  output  can  best  be  shown  by  a sample  case.  Specific  engineering  aspects  of  this 
sample  are  given  in  section  4.4  of  the  engineering  analysis  (vol.  I). 
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5.1  SAMPLE  CASE  DESCRIPTION 

This  sample  provides  an  estimation  of  community  noise  for  an  EOW  configured  aircraft. 
The  powerplant  considered  is  a typical  high  bypass  ratio  turbofan.  This  sample 
corresponds  to  the  example  illustrated  in  section  4 4 of  the  engineering  analysis  (vol.  I 
of  this  report).  Further  detail  on  the  engineering  aspects  of  this  sample  is  reported  in 
the  section  noted.  The  objective  here  is  to  illustrate  output  and  some  particulars  on  use 
of  the  computer  program 

5.1.1  INPUT 

This  sample  is  composed  of  three  cases: 

o Community  noise  estimation  for  an  EOW  aircraft  configuration  using  the  English 
unit  input  data  option 

o Same  as  above,  except  using  the  S I.  unit  input  data  option 

o Same  as  above,  except  for  an  unsuppressed  engine-under-w  ing  iEUWi  aircraft 

configuration 

The  first  two  cases  use  new  prediction  capability,  whereas  the  last  case  uses  the 
capability  originally  available  in  the  software  of  contract  NAS2-6969.  The  inpu!  data 
cards  for  this  sample  are  listed  in  section  5.2. 

A few  subtleties  need  to  be  brought  out  about  the  data  input  shown.  The  noise 
components  considered  in  the  sample  are  primary  and  secondary  jet.  turbine,  exit  fan. 
jet/edge  interaction,  inlet  fan  and  core  noise.  The  data  for  these  noise  components  are 
input  as  though  the  noise  is  produced  by  two  different  engine  wing  configurations: 

o Engine/wing  configuration  1 considers  primary  and  secondary  jet.  turbine,  exit  fan. 
and  jet/edge  interaction  noise  components 

o Engine/wing  configuration  2 considers  inlet  fan  and  core  noise  components- 

Actually,  only  one  engine  configuration  is  being  considered  although  the  data 
corresponds  in  format  as  though  there  are  two.  This  manner  of  input  is  necessary  to 
overcome  some  limitations  noted  previously: 

o Section  4.3.3  for  variable  ISW3.  Only  one  set  of  empiricisms  can  be  used  in  the 
wing  shielding  calculations  each  time  the  core  and  turbine  noise  estimation  module 
is  called. 

Section  2. 4. 1.2  of  the  engineering  analysis  (vol.  I of  this  report)  An  error  will 
result  if  an  attempt  is  made  to  estimate  the  inlet  fan  noise  diffracting  about  the 
wing’s  trailing  edge.  The  error  is  corrected  by  setting  the  variable  IWED(  1 ) of 
section  4.4  equal  to  2ero.  This  error  arises  because  of  the  engine/wing  geometry 
specified  indicates  a reflection  situation  concerning  the  inlet  position  relative  to 


i 


i 
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the  half-plane  approximating  the  trailing  edge.  No  reflection  situation  actually 
exists.  The  signal  diffracting  about  the  trailing  edge  is  in  fact  negligible  compared 
to  the  total  radiated  inlet  noise. 


After  the  input  format  for  the  data  was  selected  to  overcome  the  noted  limitations,  the 
order  of  input  for  each  noise  component  ($NOISIN  cards)  of  each  engineiwing 
configuration  is  maintained  for  all  three  cases  according  to  the  instructions  given  in 


section  4 3 


5.1.2  OUTPUT 

This  sample  produced  a large  volume  of  printed  output  because  all  output  reports  were 
selected  isee  vanab'  lOl'T  of  sec  4.1)  The  objective  here  is  to  show  the  pertinent 
format  for  each  type  of  output  report  one  can  obtain,  as  outlined  in  section  5.1  An 
example  of  the  printed  output  is  given  in  section  5.3.  In  this  example,  only  the  second 
case  in  S I units  is  illustrated  To  reduce  this  volume  of  printed  material  further, 
predictions  for  only  one  noise  component  and  one  observer  position  for  the  case  are 
shown  A complete  printout  for  the  sample  can  be  obtained  by  simply  executing  the 
computer  program  using  the  sample  input  data  shown. 
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5.2  SAMPLE  INPUT 


SHIELDING  OFCKCUT  TEST  EPR=5.  T.O.  P OWE  K ( E NGL  I SH  UMTS) 

JGOAT A AMACH=.2S,  I A T MGS  = 4 ,C PF E S = 1 A.  7,CPHUMD=  7D.0,CTENP=536.7,EPP=0., 
I0CP  = 2 , !CUT<  1 1 = 1 .2 ,3 .4.5 ,6,7 ,FLR  = 0. , 1SPTRM  = 0,  IUNIT=1, 

N08S  = 2 ,NTENG  = 2,SLDISTI 1 )=0.,2128. , ALTPC= 1000. , INSECWI l ) = 1 , 1 ,0$ 
JFWDATA  SwPTE=  10  . , SwPL  E = 2B.  .DIHfc  r = 3.  . Cl)SD  = 8.  1 .COX  10=  3.  I ,CDX20=1.5  , 
OCX  CD- 0.  . CCY1;;>  = 0. 84,CDY2D  = C.2  0.DOYU0  = .75 ,0010=  4.0,  C lANfc  = 6 .42$ 
JETIP  ♦ S),  TURBINE,  AFT  FAN, JET  ECGE  I NT  f P ACT  I ON  COMPONENT  S 
$ NO  I SIN  MTYPt  = A,NENG  = 2, I TYPE =2 , N JE T I = 3 . MCGCE 1 =2 , 

AP1=6.65,VP1  = 1370.,,.PI  = 263.,AS2  = 20.G4,VS2  = 919.9,WS2=1312., 
INS.HL0  = I ,LIAMT1  = 3.96$ 

$ SHL DAT  OSL 1=5.36,CSL2=3.32,DIAMT2=&.42$ 

JNOISIN  I T yp  E = 3 , I Sv*3  = 3 , IC3=0.BN3=1  10. , CS 3 = .6 6 , PMF 3= 26 3 . , 

SS3  = 3404.  , TU2=  1478., VTR  3=726.9$ 

$ SHL OAT  EMJ=0.A6,TSTS(’  = 2.56,  I wED  (1)  = 1 . 1 , l $ 

JNO  I SIN  I TYPE  = 5,NSTGA5  = l ,FPRA5( I) =1. 55,NQA5( 1 ) = A6,PN1A5=3A0A., 

R SS45I  II  =?CO.,PT  Sr5=  1 . 34  ,7.ce  A5<  1 ) = 20.04,6PR5  = 5.  ,NI  5 = 2$ 
t SHL DAT  £MJ =C. 82, TSTSC= 1 .05, I WEDI  1 > = l , 1, 1 * 

INC  I S IN  ITYPE= 13  .ECGVAP ( l ) = 3.l  ,.62,6.42,1.05,C.,5.7,.7854,.25,5.,O.J 
S E V,  U A T A SWT  TE=  in.  , SWPL  f = 2H.  , D I HE  rj  = x.  , 0 0,  S 0 = 0 - 1 .COX  1 D = . . I ,00X20=  1.5  , 

C 0x00  = 0.  , Cf  Y 10  = 0.84,  DO  Y2D  = C.  2 8,  !>OYOD  = . 75  , DDL  0 = A . 0 , 0 I ANF  = 6.  A2  $ 
INLET  FAN  ANC  CC^E  NOISE  COMPONENTS 
JNO I SIN  NTYP£  = 2,NEN0  = 2, I TYPE  = 4,  N S TC-4  1>=  1 ,FPC/,5 ( 1 1 = 1 . 55, NBA  5 ( 1)=A6, 
RMA5--3ACA.  , c SSA6U  ) =300.  ,t  TSA5=  l.  3A,  CFPPA  = 1.29, 

L INA=l,CFA  = l 116. , E DH A = 3 • 6 , t L UHA= 1 . 5 , F MA  = -0 . A , I OP A= 2 , 

ILAYA  = 2,Nv>LA  = l, GIAMBI  1)  = 7.69,  I NSHLD=  1 i 
JSHLCAT  IWCC!1)=C,1, 1, IMJSP=OS 

JNO  IS  IN  I T YFE  = 3,  IS, ,3=2  , IMF  3 = 2 9.2, EK3  = 2. ,JB3=1 ,PP3= 13.94, TT 3 = 2 795  .$ 
JSHLCAT  EM J =0.46  , T ST SJ  = ? .55, I wEl  (11=1,1,1$ 

SH I ELDING  CHECKOUT  TEST  PPP=5,  T.O.  PCwFR(MKS  UNITS) 

J GO AT  A CPE t S = l. C.CRHU MU =70. ,CTF  MP  = 2 96.2 , 1 UN  I T=0, 
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6.0  MACHINE  REQUIREMENTS 


This  program  is  designed  to  operate  on  a CDC6600  scientific  computer.  Approximately 
141. IK  octal  words  of  storage  are  required  for  operation.  Data  input  is  thorugh  cards, 
tape  or  disk  card  image.  Output  is  to  a line  printer.  In  addition  to  the  standard 
input/output  disk  files  TAPE5/TAPE6,  the  following  disk  files  are  used  as  scratch 
storage  (primarily  report  files)  TAPE8,  9,  10,  11,  12,  13,  and  20.  TAPE20  is  the  only 
scratch  file  that  the  user  may  need  to  concern  himself  with.  This  file  is  used  for  output 
of  acoustic  data  for  noise  contour  estimation  (see  NASA  CR-114649,  pp.  38,  39.  175  to 
178  for  further  detail). 


7.0  OPERATING  SYSTEM 


The  program  has  been  checked  out  on  the  CDC6600  using  the  FTN  4.5  compiler  under 
the  KRONOS  2.1  operating  system.  The  program  is  written  in  FORTRAN  IV  language 
to  be  relatively  machine  independent.  The  majority  of  the  FORTRAN  code  complies 
with  ANSI  standards,  except  for  the  use  of  NAMELIST  and  OVERLAY  features 
designed  into  the  program  for  economical  and  easy  use. 


8.0  RESOURCE  ESTIMATES 


The  central  processor  unit  (CPU)  time  required  to  process  a job  depends  upon  the 

program  options  used.  The  major  factors  influencing  the  time  per  case  are: 

1.  The  number  of  configurations  and  the  number  of  noise  components  in  each 
configuration 

2.  1/1  or  1/3  octave  bands  for  predicted  noise  spectra  (the  1/3  octave  band  option  uses 
approximately  twice  as  much  CP  time) 

3.  The  number  of  sideline  observer  positions 

4.  The  number  of  noise  components  for  which  shielding  is  included,  and  the  number 
of  wing  edges 

5.  Lining  attenuation  and  configuration  corrections 

6.  Optional  output  reports 

The  execution  time  in  CPU  seconds  per  case  is  given  approximately  by 

N1  N2  (Cl  + C2  N3  + N4(C3  N3  + C4)) 


where 

N1  = number  of  noise  components 
N2  = number  of  observer  sideline  positions 
N3  = number  of  wing  edges 
N4  = number  of  frequency  bands 

(Cl,  C2,  C3,  C4)  = (1.0,  0.3,  0.02,  0.04) 

This  formula  is  based  on  computer  runs  with  all  print  options  on  the  shielding  included 
for  each  noise  component  (for  jet  noise  shielding  assume  N3  = 3,  for  jet  edge  N3  0i. 
For  lining  attenuation  and  configuration  corrections,  add  approximately  ](y?  to  the 
execution  time. 


9.0  DIAGNOSTICS 


i 
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i 
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The  following  is  a list  of  the  diagnostic  messages  which  are  printed  when  various  error 
conditions  are  detected  by  the  program. 


1.  TOO  Many  entries  in  altitude  vs  temperature  table. 

2.  TOO  MANY  ENTRIES  IN  ALTITUDE  VS  PRESSURE  TABLE.  MAXIMUM 
; Lf.OWED  IS  FIFTY.  IE,'  ATMOSPHERE  IS  ASSUMED. 

3.  TOO  MANY  ENTRIES  IN  ALTITUDE  VS  RELATIE  HUMIDITY  TABLE. 
MAXIMUM  ALLOWED  IS  FIFTY.  ISA  ATMOSPHERE  IS  ASSUMED. 

4.  ALTITUDE  VS  TEMPERATURE  TABLE  IS  UNDEFINED.  MUST  HAVE  AT 
LEAST  TWO  ENTRIES.  ISA  ATMOSPHERE  IS  ASSUMED. 

5.  ALTITUDE  VS  PRESSURE  TABLE  IS  UNDEFINED.  MUST  HAVE  AT  LEAST 
TWO  ENTRIES.  ISA  ATMOSPHERE  IS  ASSUMED. 

6.  ALTITUDE  VS  RELATIVE  HUMIDITY  TABLE  IS  UNDEFINED.  MUST  HAVE 
AT  LEAST  TWO  ENTRIES.  ISA  ATMOSPHERE  IS  ASSUMED. 

7.  EFFECTIVE  TIP  MACH  NUMBER  OUT  OF  RANGE  <GT  0.93  OR  LT  O.i  OR  BAD 
INPUTS. 

3.  TOO  MANY  TARGET  FREQUENCIES  SPECIFIED  FOR  LINING.  ONLY  FIRST 
TEN  ARE  USED 

9.  TOO  MANY  WALLS  SPECIF  ED  IN  FAN  LINING.  ONLY  FIRST  TEN  ARE 
USED. 

10  NO  WALLS  HAVE  BEEN  DEFINED  FOR  FAN  LINING. 

11.  ERROR  WR'TING  RANDOM  FILE  JOB  ABORT. 

12.  ERROR  READING  RANDOM  FILE  JOB  ABORT. 

13.  *** FATAL  ERROR*** BAD  GEOMETRY. 

14.  DO  CORE  SEPARATE  FROM  TURBINE  WHEN  SHIELDING  INCL  USED 
TURBINE  ASSUMED. 
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10.0  CONTROL  CARDS 


A computer  magnetic  tape  (with  a backup  tape)  is  furnished  for  the  complete 
NAS2-6969  contract  software  revised  in  accordance  with  the  requirements  of  contract 
DOT-FA74WA-3497.  The  tape  contains  the  files  listed  in  the  sequence  below  created  by 
the  text  program  editor  UPDATE: 

SOURCE  . . File  of  source  deck  for  regeneration  of  NEWPL  and  COMPILE  files 
NEWPL  . . .New  program  library  file 

COMPILE  . . .Card  image  file  for  compilation  using  the  FTN  4.5  compiler 

Any  of  these  files  on  the  tape  can  be  used  to  create  a file  of  relocatable  binaries  to  be 
executed  on  the  CDC6600  computer. 

The  software  on  the  tape  actually  corresponds  to  three  independent  programs;  i.e., 

1.  A community  noise  source  estimation  program  revised  to  include  airframe 
shielding  attenuation  estimation  capability 

2.  A postprocessor  program  to  interface  the  output  of  item  1 with  a noise  contour 
estimation  program 

3.  A noise  contour  estimation  program 

The  first  program  is  the  only  one  which  was  significantly  modified  for  the  DOT  FAA 
contract.  The  other  two  programs  are  just  °DC6600  versions  of  the  original  IBM360 
computer  versions  delivered  to  NASA-Ames  in  July  of  1973.  Use  of  the  last  two 
programs  follows  that  described  in  NASA  reports  CR114649  and  CRl  14650  except  for 
the  following: 

4.  The  $ notation  for  NAMELIST  input,  on  the  CDC6600  computer  replaces  the  & 
notation  used  for  the  IBM360  computer. 

5.  The  control  cards  are  radically  different  for  the  CDC6600  computer  for  all  the 
programs  mentioned  above. 

There  are  many  ways  to  get  the  software  code  off  of  the  tape(s)  provided  We  propose  the 
following  methods  which  result  in  economical  and  easy  use  for  a large  number  of  user's. 
It  is  assumed  that  the  CDC6600  computer  system  has  indirect-access  permanent  file 
disk  storage  capability.  The  storage  space  required  for  all  the  programs  is  about  800 
sectors  (50K  words  or  3M  bits) 

The  control  cards  sequences  used  below  are  unique  to  the  Boeing  Enhanced  KRONOS 
2.1  computer  system,  but  these  can  be  used  at  other  CDC6600  computer  installations 
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with  minor  modifications.  These  control  card  sequences  require  certain  information  to 
be  provided  by  the  user  as  indicated  in  table  A-2. 

The  first  step  that  must  be  accomplished  is  to  get  the  programs  off  of  the  tape,  compiled, 
and  retained  on  disk  storage  as  relocatable  binary  programs.  The  control  card  sequence 
listed  in  table  A-3  instructs  the  computer  to  perform  this  task.  Details  involved  are 
noted  in  the  table  in  the  description  column. 

Once  the  relocatable  binary  programs  are  retained  on  disk,  the  control  card  sequences 
given  in  tables  A-4,  A-5,  and  A-6  can  be  used  to  load  and  execute  the  programs. 
Table  A-4  applies  for  the  noise  source  estimation  program.  Table  A-5  applies  for  the 
postprocessor  program  to  interface  acoustic  data  output  of  a previous  noise  source 
estimation  job  with  the  noise  contour  estimation  program.  Table  A-6  applies  for  the  use 
of  the  noise  contour  estimation  program. 


Abbreviation 


yy 

Userno 


Password 
Name,  etc 


--  eoi  -- 


Table  A-2.  — User  Supplied  Information 

Definition 

Unique  job  name 
^ob  priority  selected  by  user 
User's  account  number 
User's  password 

Users  name/telephone  number/mail  stop/ 
organization  number 
User's  mail  stop 

Tape  identification  number  for  source 
program  library  tape 

Central  processor  time  estimate  in  seconds 
Estimate  of  the  maximum  number  of  cards 
punched.  Usually  this  is  less  than  500 
cards  per  job  step. 

Permanent  file  name  for  Noise  Source 
Estimation  (NSE)  overlay  program  binaries 
Permanent  file  name  for  Post-Processor  (PP) 
program  binaries 

Permanent  file  name  for  Noise  Contour 
Estimation  (NCE)  main  program  binary 
Permanent  file  name  for  Alternate  LIBrary  (ALIB) 
containing  subroutine  binaries  used  by  NSE  and 
NCE  programs 

Permanent  file  name  to  save  computed  noise 
contour  data  for  subsequent  run  to  produce 
plots 

End  of  record  card;  i.e.,  a card  having  a multiple 
789  punch  in  column  one  with  any  informative  alpha- 
numeric characters  punched  on  the  rest  of  the  card. 
End  of  information  card;  i.e.,  a card  having  a 
multiple  6789  punch  in  column  one  with  any  infor- 
mative alphanumeric  characters  punched  on  the  rest 
of  the  card. 
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Table  A-3.— Computer  Job  to  Save  Relocatable  Binaries 
on  Disk  From  Program  Library  Tape 


Control  Cards 


Description 


Jobn  ,CM120000,T200,Pyy. 

ACC0UNT,Userno, Password.  Name,  etc. 
RFL, 20000. 

REQUEST, PL,  VSN=66xxxx,F=I,LB=KU,P0=AR. 


Job  card 
Account  card 


REWIND, PL, INPUT. 

SKIPF  ,PL  ,1 . 

C0PYNF.PL ,0LDPL  ,1 . 

RETURN.PL. 

C0PYSBF, INPUT, 0UTPUT. 

REWIND, INPUT. 

SKIPR, INPUT,  1. 

RFL, 40000. 

UPDATE ,Q,C. 

RFL, 120000. 

FTN , I=C0MPILE ,T ,PL=77777B. 

RFL, 20000. 

REWIND.LG0. 

C0PYBR,LG0,A,15. 

C0PYBR ,LG0,B ,2. 

C0PYBR  ,LG0 ,C  ,1 . 

C0PYBF,LG0,D. 

REWIND, D. 

, . LIBGEN ,F=D,P=ALIB,N=ALIB. 

; i SAVE ,A=pfnl/CT=S ,M=R. 

; | SAVE ,B=pfn2/CT=S ,M=R.  I 

SAVE ,C=pfn3/CT=S ,M=R. 

SAVE  ,ALIB=pfn4/CT=S ,M=R.  __ 

--  eor  -- 

f *C0MPI LE  PHBNAC . ZERO 

S'  ----- 

r t 

l 
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Gets  program  library  tape  and  assigns 
file  name  PL  to  the  tape 


Copies  program  library  tape  to  0LDPL  file 
Copies  input  deck  to  0UTPUT  file 


Writes  program  "decks"  on  C0MPILE  file 

Compiles  program  to  produce  relocatable 
binaries  on  LG0  file 

Copies  NSE  overlay  programs  to  file  A 

Copies  PP  program  to  file  B 

Copies  NCE  main  program  to  file  C 

Copies  subroutines  for  NSE  and  NCE 
programs  to  file  D 

Generates  an  alternate  library  on  file 
ALIB  for  loader  from  routines  on  file  D 

Saves  files  A,  B,  C,  ALIB  on  disk 
for  later  use 

Instructs  UPDATE  editor  to  write  "decks" 
PHBNAC  through  ZER0,  to  the  C0MPILE  file 


Table  A-4.— Computer  Job  to  Execute  Noise  Source  Estimation  Program 


Control  Cards  

Jobn ,01141100, Txxx ,Pyy. 

ACCOUNT, Userno, Password.  Name,  etc. 
RFL, 20000. 

GET,NSE=pfnl/UN=Userno. 

GET ,ALIB=pfn4/UN=Userno. 

OFFLINE.  PUNCH  EST  = ZZZ  CARDS3 
REWIND, INPUT. 

C0PYSBF, INPUT .OUTPUT. 

REWIND, INPUT. 

SKIPR, INPUT, 1. 

RFL, 141100. 

L0ADXEQ,F=NSE,U=ALIB. 

RFL ,20000. 

C0PYCF , TAPE20 , PUNCH ,1,1,80. 


RFL, 141100. 

EXEC. 

RFL, 20000. 

COP YCF .TAPE20 , PUNCH ,1,1,80. 

--  eor  -- 

Data  deck  for  1st  job  step 
--  eor  -- 

Oata  decks  for  additional  job  steps 
--  eoi  -- 


Description 
Job  card 
Account  card 


Gets  NSE  program  off  disk 

Gets  alternate  subroutine  library  off  disk 

Need  only  if  cards  are  to  be  punched 

Writes  input  deck  to  OUTPUT  file 


Loads  and  executes  NSE  program 


Need  only  if  acoustic  data  is  punched 
for  noise  contour  estimation 


Repeat  these  cards  to  execute  NSE 
program  for  each  additional  job  step; 
i.e.,  you  have  more  than  one  data  deck 


Repeat0  the  data  decks  for  each 
additional  job  step 


a)  Cards  punched  are  acoustic  data  (noise  level,  engine  performance 
parameter,  elevation  angle,  log^Q  of  range  at  CPA)  which  characterize 
an  airplane  configuration.  See  NASA  CR114649,  pp.  38,  39  and  175-178. 
Format  for  punch  data  on  the  cards  is  (1PE12.3,  3E12.3). 

b)  Not  needed  if  only  one  data  deck  is  used. 


Table  A-5.— Computer  Job  to  Execute  Postprocessor  Program 


Control  Cards 


Jobn ,CM60000,T100,Pyy.  v 
ACC0UNT,Usemo, Password.  Name,  etc. 
RFL, 20000. 

GET,PP=pfn2/UN=Usemo. 

0FFLINE.  PUNCH  EST  = zzz  CARDS3 
REWIND, INPUT. 

C0PYSBF, INPUT, 0UTPUT. 

REWIND, INPUT. 

SKIPR, INPUT,  1. 

...  —> 

RFL, 33000. 

REWIND .TAPE20. 

C0PYCR, INPUT, TAPE20.1 ,1 ,80. 

PP. 

RFL,  60000. 

FTN , I=TAPE22 ,T ,B=PUNCHB . 

--  eor  -- 

Data  deck*5  for  1st  job  step 
--  eor  -- 

Data  decksb  for  additional  job  steps 

— -1 

--  eoi  -- 


Description 


Job  card 
Account  card 


Gets  post  processor  program  off  disk 


Copies  input  deck  to  0UTPUT  file 


Produces  acoustic  data  subroutine  to 
interface  with  NCE  program,  compiles 
the  subroutine,  and  punches  relocatable 
binary  deck.  Repeat°these  cards  for 
each  additional  job  step  to  characterize 
more  than  one  airplane  configuration. 


Data  deck  of  acoustic  data  for 
NCE  program 

Include0  additional  data  decks  for 
characterizing  the  noise  of  more 
than  one  airplane  configuratioh. 


Cards  punched  are  binary  acoustic  data  subroutines  to  be  used 
by  the  noise  contour  estimation  program. 

Input  data  cards  are  acoustic  data  (noise  level,  engine  performance 
parameter,  elevation  angle,  log1Q  of  range  at  CPA)  which  characterize 
an  airplane  configuration.  See  NASA  CR114649,  pp.  38,  39  and  175-178. 
Format  for  data  on  the  cards  is  (1PE12.3,  3E12.3). 

Not  needed  if  only  one  data  deck  is  used. 


Table  A-6.~Computer  Job  to  Execute  Noise  Contour  Estimation  Program 


Control  Cards 


Jobn .CM70000 ,T200,Pyy. 
ACC0UNT,Userno, Password.  Name,  etc. 
RFL, 20000. 

GET,C=pfn3/UN=Usemo. 

GET,ALIB=pfn4/UN=U$erno. 

GET , TAPE2=pfn5/UN=Userno . 
C0PYBR.C.NCE.1. 

C0PYBR,  INPUT  ,NCE,1. 

C0PYSBF, INPUT .0UTPUT. 

REWIND, INPUT. 

SKIPR,INPUT,2. 

UNBL0CK.TAPE99. 

L0ADXEQ , F=NCE , U=ALIB .CALCMPF . 
REPLACE. TAPE2=pfn5. 

PL0TFIL .CALC0MP ,TAPE99,0. 

C0MMENT.  PLAIN  WHITE  PAPER 
C0MMENT . BLACK  WET  INK 
C0MMENT.  RIGHT  EDGE  START 
C0MMENT.  MAIL  T0  M/S  xx-xx 


Description 


Job  card 
Account  card 

Gets  NCE  main  program  off  disk. 

Gets  alternate  subroutine  library  off  disk 

Getsc  previously  saved  noise  contour 
data  off  disk. 

Includes  acoustic  data  routine  with 
NCE  program. 

Copies  input  data  to  0UTPUT  file. 


Loads  and  executes  NCE  program  using 
alternate  libraries  ALIB  and  CALCMPF. 
Saves0  noise  contour  data  for  later  use. 

Disposes  plot  file  TAPE99  for  offline 
CALC0MP  plotting. 


Binary  deck  of  acoustic  data  subroutine3 
— eor  — 

Input  data  cards^  to  describe  flight 
paths,  noise  contours  desired,  and/or 
plotting  options 
--  eoi  -- 


One  of  the  routines  produced  by  the  post-processor  program  (Table  5). 

See  NASA  CR114649,  pp.  175-182  and  NASA  CR114650,  pp.  66-108  for  input 
data  description.  Note  that  a $ replaces  the  & notation  given  in  the 
second  reference  for  NAMELIST  input. 

Required  only  if  computed  results  are  being  checked  before  making  plots 
(see  NASA  CR114650,  pp,  60  and  7Q  for  discussion)  with  a second  job  submittal. 


f 

f. 


11.0  SOURCE  LISTING 


The  source  listing  for  the  programs  on  magnetic  tape  is  provided  in  volume  III  of  this 
report. 
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APPENDIX  B 

PRELIMINARY  TEST  PLAN  FOR  ASSESSMENT  OF  FORWARD 
VELOCITY  EFFECTS  ON  THE  SHIELDING  EFFECTIVENESS 
OF  WING-  AND  BODY-TYPE  ACOUSTIC  BARRIERS 


J.  M.  Campbell 


1.0  SUMMARY 


A test  plan  providing  a logical  follow-on  to  theory  and  analysis  reported  in  volume  I is 
presented.  Tests  executed  in  a wind  tunnel  using  simple  geometric  shapes  in  a variety 
of  configurations  and  for  a range  of  tunnel  flow  velocities  are  recommended  The 
proposed  matrix  of  configurations/tunnel  flow  velocities  is  aimed  at  isolating  and 
assessing: 


o Effect  of  engine  exhaust  flow  on  shielding 

o Effect  of  ambient  flow  velocity  on  the  intrinsic  shielding  processes 
o Effect  of  wing  wakes  on  wing  shielding  effectiveness 

Special  pulse  test  techniques  and  associated  instrumentation  are  proposed  to  eliminate 
problems  of  background  noise  and  reverberation  which  have  typically  dilluted  the 
usefulness  of  previous  small  static  chamber  and  wind  tunnel  data. 


2.0  INTRODUCTION 


The  mathematical  models  developed  during  this  contract  for  predicting  shielding  effects 
are  based  on  data  gathered  from  a variety  of  sources.  These  data  were  derived  from 
small  scale  static  tests  and  full  scale  static  and  flyby  tests.  In  some  cases,  these  tests 
were  only  peripherally  directed  toward  shielding  effects.  The  data  in  many  cases  are 
contaminated  by  acoustic  energy  reflected  from  the  ground  or  other  surfaces  and  noise 
sources  additional  to  the  shielded  source.  The  reflection  of  energy  from  the  ground  plane 
or  any  other  flat  surface  is  not  uniform  at  all  frequencies.  Phenomena  such  as  wind, 
temperature,  humidity,  and  their  gradients  can  change  the  amplitude  and  frequency 
response  of  the  reflection  effects.  All  full  scale  data  are  subject  to  multiple  source 
interference. 

New  data  acquisition  methods  are  needed.  Data  should  be  taken  in  a new  way.  Boeing 
Commercial  Airplane  Company,  while  improving  its  products,  has  developed  and 
reduced  to  practice  this  needed  new  way.  The  test  plan  presented  here  makes  full  use  of 
the  most  advanced  techniques  in  acoustic  measurement  to  obtain  noise  shielding  data  in 
a wind  tunnel  uncontaminated  by  reflections  and  stray  noise  sources. 


3.0  TEST  OBJECTIVES 
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The  test  program  described  herein  is  designed  and  sequenced  to  answer  three  general 
questions: 

• What  effect  does  airflow  have  on  the  "shielding”  provided  by  aerodynamic  lifting  or 
control  surfaces? 


For  some  specific  geometries,  the  aircraft  body  provides  "shielding.”  How  is  this 
affected  by  airflow? 


• Does  a wake  provide  significant  "shielding”  in  a flow  field? 

The  term  "shielding”  includes  reflection,  absorption,  refraction,  and  diffraction  in 
differing  proportions  depending  on  the  situation.  By  programming  the  tests  properly, 
these  proportions  can  be  revealed.  Thus,  the  quality  as  well  as  the  quantity  of  the 
"shielding”  may  be  determined.  While  flight  effects  cannot  be  measured  directly  in  a 
wind  tunnel,  this  test  program  is  a necessary  step  in  interpreting  measured  flight  data. 


1 1 
[ 


[ i 

y 

I i 


I < 


fc 


I 


4.0  TEST  APPROACH 


The  two  principal  problems  of  wind  tunnel  acoustic  measurements  are:  high  noise  levels 
in  the  tunnel  (see  fig.  B-l)  and  the  hard  walls  of  the  tunnel  provide  reflecting  surfaces 
which  cause  a reverberent  level  to  build  up.  The  solutions  to  these  problems  are 
interrelated  so  that  the  "best”  solution  for  one  is  the  "worst”  for  the  other. 
Nevertheless,  a viable  compromise  is  available. 

The  effect  we  wish  to  measure  is  shielding.  Shielding  can  be  defined  as  excess  sound 
attenuation  between  two  points,  excess  in  that  the  attenuation  is  more  than  that  which 
is  expected  from  inverse  square  losses  and  atmospheric  absorption.  Since  attenuation  is 
the  measurement  desired,  the  absolute  magnitude  is  not  an  important  parameter-only 
the  change  is. 

The  recommended  test  hardware  arrangement  in  the  wind  tunnel  is  shown  in 
figure  B-2.  A sound  source  is  placed  in  position  on  one  side  of  a shielding  configuration, 
and  a microphone  is  placed  close  to  the  source.  This  microphone  is  called  the  reference 
microphone.  A second  microphone,  called  the  receiver  microphone,  is  mounted  at  a 
position  of  interest.  If  the  signal  from  the  reference  microphone  is  compared  to  the 
receiver  signal,  the  change  will  represent  the  total  attenuation  between  the  two  points. 

If  the  signal  from  the  sound  source  is  pulsed,  and  the  signals  from  the  two  microphones 
are  examined  before  any  significant  energy  is  reflected  from  the  wind  tunnel  walls,  the 
data  will  be  free  from  the  effects  of  reverberation.  Using  the  Boeing  transonic  wind 
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THIRD-OCTAVE  BAND  LEVEL  IN  dB  RE  20  (iN/M2 


SOUND  SOURCE 


Figure  B 2.  General  Test  Hardware  Arrangement  in  Wind  Tunnel 


tunnel  (BTWT)  as  an  example,  the  source  and  reference  microphone  would  be  placed 
0.60  m above  the  test  article,  while  the  receiver  would  be  mounted  0.60  m below 
(fig.  B-2). 

The  path  length  around  the  test  article  is  1.7  m.  The  transit  time  would  be  about  5 ms. 
The  shortest  reflection  path  would  be  about  1.2  m longer  than  the  most  direct  path; 
therefore,  the  longest  time  available  to  examine  the  two  signals  is  about  4 ms. 

The  effect  of  high  noise  levels  in  the  wind  tunnel  (fig.  B-l)  would  be  alleviated  by  using 
a high  intensity  sound  source  which  projects  a single  tone.  The  choice  of  the  frequency 
of  the  tone  is  dependent  on  the  wind  tunnel  used.  Referring  to  figure  B-l.  the  sound 
level  in  the  BTWT  in  the  1/3  octave  band  centered  about  16  kHz  is  73  dB.  Consider  the 
sound  source  is  set  on  this  frequency  and  produces  a level  of  120  dB  at  1 m.  In 
traversing  the  1.7-m  path  to  the  receiving  microphone,  the  normal  attenuation  would  be 
5 dB,  resulting  in  a level  of  115  dB.  We  would  expect  shielding  effects  up  to  40  dB; 
hence,  the  level  of  the  signal  at  the  receiving  microphone  would  be  as  low  as  75  dB,  only 
2 dB  above  the  ambient.  A solution  to  the  poor  signal-to-noise  ratio  is  to  limit  the 
receiver  bandwidth.  Referring  to  the  previous  paragraph,  we  find  that  the  longest 
reverberation  free  time  is  4 ms.  With  this  length  of  time,  the  narrowest  bandwidth 
achievable  is  1 (4xl0'3)  or  250  Hz.  Considering  that  the  13  octave  band  about  the 
16-kHz  center  frequency  is  about  3750  Hz  wide,  the  ambient  noise  reduction  would  be 
10  log,„i3750/250)  or  -12  dB,  giving  us  a signal-to-noise  ratio  of  14  dB.  This  would 
introduce  a systematic  error  of  only  0.3  dB  for  a shielding  value  of  40  dB.  For  shielding 
values  of  34  dB  or  less,  the  error  becomes  less  than  0.1  dB.  If  the  source  frequency  is 
raised  to  a higher  frequency,  the  percentage  bandwidth  reduction  could  be  greater,  but 
the  atmospheric  absorption  also  becomes  greater,  thus,  cancelling  the  improvement 
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5.0  TEST  PROCEDURE 


Each  configuration  will  be  tested  in  the  wind  tunnel  for  wind-off  and  three  wind-on 
conditions.  A configuration  is  defined  as  a particular  combination  of  test  article,  sound 
source  placement,  and  receiver  placement.  The  neutral  airfoil  qualifies  as  the  most 
tested  article  (table  B-l).  It  will  be  used  to  establish  the  basic  wind-off/wind-on 
parametes  for  both  leading-  and  trailing-edge  shielding.  It  is  also  used  to  gain  frequency 
dependent  data.  The  wing  is  to  be  tested  not  only  by  itself,  but  also  in  combination  with 
simulated  engine  flow.  The  round  nozzle  and  the  D-shaped  nozzle  will  be  tested  alone 
and  with  the  wing.  The  sound  source  will  be  placed  inside  the  nozzle  for  these  tests.  A 
cylinder  will  be  tested  alone  to  simulate  a fuselage.  In  addition  to  tests  with  a test 
article  in  the  tunnel,  a series  of  ambient  noise  tests  will  be  made  and  an  additional 
series  with  the  sound  source  activated.  la  tat*l,  495  configurations  will  be  tested.  This 
corresponds  to  1980  data  points. 
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6.0  TEST  CONFIGURATIONS 


It  is  planned  to  test  in  the  wind  tunnel  eight  configurations: 

1.  A thin  neutral  airfoil 

2.  A wing  section 

3.  A round  convergent  nozzle 

4.  A D-shaped  nozzle 

5.  A round  convergent  nozzle  below  a wing  section 

6.  A round  convergent  nozzle  above  a wing  section  (unattached  flow) 

7.  A D-shaped  nozzle  above  a wing  section  (attached  flow) 

8.  A cylinder 

Sketches  of  these  configurations  will  be  found  in  figure  B-3. 
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7.0  ACOUSTIC  INSTRUMENTATION  AND  DATA  REDUCTION 

The  acoustic  equipment  consists  of  two  groupings  (fig.  B-4).  One  grouping  is  locate.d  in 
or  adjacent  to  the  wind  tu/inel;  the  second  grouping  is  located  in  building  2-80  at  The 
Boeing  Company  Plant  II  complex.  Collectively,  the  second  set  of  equipment  is  called 
the  Acoustic  Data  Processing  Facility  (AD P-1).  The  first  set  of  instrumentation  consists 
of  the  sound  source  (including  its  modulator),  the  two  microphones,  the 
signal-conditioning  equipment,  and  monitoring  equipment.  The  ADP-1  is  an  on-line 
real-time  data  acquisition  and  reduction  system.  The  heart  of  the  ADP-1  is  the 
minicomputer  manufactured  by  Prime,  Inc.  This  computer  is  a multilingual 
time-sharing  system.  It  has  a basic  memory  of  64,000  words,  and  the  memory  can  be 
expanded  to  256,000  words  in  the  future.  The  block  diagram  in  figure  B-5  indicates  the 
capabilities  of  the  system.  For  these  tests  only,  a portion  of  this  capability  will  be 
employed. 

The  equipment  would  be  used  in  the  following  way.  Upon  command  of  the  operator,  the 
ADP-1  will  provide  a 10-ms  pulse  to  the  sound  source  modulator  (fig.  B-4).  This  causes 
the  sound  source  to  generate  a tone  burst  at  a selected  frequency.  The  reference 
microphone  is  placed  so  as  to  pick  up  the  acoustic  signal  almost  immediately.  The 
receiver  microphone  captures  the  signal  at  some  later  time.  Both  signals  are  digitized 
and  are  processed  through  a "digital  filter.”  The  amplitude  of  each  pulse  as  well  as  the 
ratio  between  them  is  determined  and  recorded.  As  soon  as  this  process  is  complete, 
ADP-1  issues  another  pulse  and  the  same  process  is  executed.  Upon  completion  of  30 
cycles  of  this  data  gathering,  ADP-1  calculates  the  mean  and  standard  deviation  of  the 
reference  signal,  the  receiver  signal,  and  the  ratio  of  the  individual  pulses. 
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Figure  B- 


Figure  B-5. -Schematic  of  Acoustic  Data  Processing  System  ( ADP-1 ) 


8.0  DATA  ANALYSIS 


The  data  output  will  consist  of  the  mean  and  standard  deviation  of  the  reference  and 
receiver  microphones  output  and  the  mean  and  standard  deviation  of  the  ratio  of  the 
individual  pulse  pairs.  By  computing  the  ratio  of  the  means  of  the  reference  and 
receiver  signals  and  comparing  this  with  the  mean  of  the  ratio,  it  is  possible  to 
determine  the  stationarity  of  the  process.  The  delay  time  between  the  two  signals  will 
show  the  path  length  and  also  the  normal  attenuation.  This  information  can  be  used  to 
examine  refraction  effects.  The  standard  deviation  of  the  reference,  as  compared  to  that 
of  the  receiver,  gives  an  additional  measure  of  the  randomness  of  the  process.  The 
target  of  the  analysis  is  comparison  between  the  data  and  the  mathematical  models 
developed  in  this  contract.  Thus,  the  mean  of  the  ratio  will  be  plotted  against  the 
significant  variables  such  as,  tunnel  wind  speed,  angle  of  source  relative  to  trailing 
edge,  angle  of  receiver  relative  to  trailing  edge,  wake  intensity  and  thickness,  etc. 


9.0  SCHEDULE  AND  COST  ESTIMATE 


A schedule  and  manpower  estimate  for  this  proposed  program  is  presented  in  table  B-2. 
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Schedule  and  Cost  Estimates 


APPENDIX  C 

UNIFORM  FLOW  EFFECTS  ON  HALF-PLANE  DIFFRACTION 

D.  G.  Dunn 


ro  = (-r0  cos  60,  r0  sin  60,  0) 

= (-x0.  y0.  0) 

T = (r  cos  0,  -r  sin  0,  0) 

K = (K  cos  0,  -K  sin  0,  0) 

M0  Parallel  to  half-plane  and 
perpendicular  to  edge 


Consider  the  case  illustrated  in  this  sketch  without  flow  where  sound  from  a source,  S, 
is  incident  upon  a semi-infinite  half-plane.  The  incident  field  for  an  omnidirectional 
compact  source  is  given  by: 


= exp[i  K • (r- r0)]/[K  * (r- r0)] 

= A(K  r)  exp(i  K rQ  cos(0  - 0Q)1 

where  A(K  r)  = exp(i  K r)/|(K  r)  [l  + (-p-)cos(0  -0O)]} 
= exp  (i  K r)/(K  r)  in  the  far-field 


(C-l) 


It  can  be  shown  by  the  application  of  the  reciprocity  theorem  illustrated  in  figure  C-l  to 
the  half-plane  that  the  Green’s  function  in  the  far  field  is  equivalent  to  the  problem  of  a 
plane-wave  incident  upon  the  half-plane  and  the  determination  of  the  field  close  to  the 
barrier.  This  latter  problem  has  been  solved  by  Sommerfelt.1  Since  we  are  interested  in 
a solution  in  the  far-field  relative  to  free-field  radiation,  the  reciprocity  theorem 
permits  us  to  drop  the  A(Kr)  factor  and  consider  the  radiation  field  as  a plane-wave 
problem;  i.e., 


% 

*, 


<t> j = expti  K rQ  cos(0  -0O)] 

- exp[i  K (-  x0  cos  6 + yQ  sin  0)1 


(C-2) 


i 
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A.  Sommerfelt,  "Optics,  Lectures  on  Theoretical  Physics,”  Vol.  IV,  Academic  Press,  Inc.,  1954. 
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The  total  field  being  sought  is 


0t  - 0j  + <p 


where  <f>  is  the  diffracted  or  scattered  field. 


The  boundary  condition  that  impacts  the  total  radiated  field  is 


■ = 0 for  y = 0,  x < 0 


The  scattered  field  0 must  satisfy  the  wave  equation 


V2  <*,  + K-  0 = 0 


As  noted,  this  problem  has  been  solved  previously.  The  solution  is 


< =G(T|70) 


= Fla  b)  + Fla  c)  expl-  i a"d) 


where  F(x)  = 


a = 2 K r„ 


r°° 

:/  exp(iz“) 


b = sin  (0  - 0o)/2 

c = sin  ( 0 + 0o)/2 


d = sin  S sin  0o 

Consider  the  case  with  flow;  i.e.,  |M0j 
and  are  to  be  replaced  by 


= 0.  Equations  (C-2)  and  (C-5)  no  longer  apply 


0j  = cxpf  i K I-  x0  cos  6 + y0  sin  0 )/( I - MQ  cos  0 )) 

+ K“  (l  +-jr  M0  • v)*-  0 = 0 or 

+^+^+i2MoKff+K-0  = O 

r)x “ ' t)y~  dt- 
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This  new  problem  is  solved  by  use  of  a modified  Lorentz  transformation  which  converts 
the  governing  equations  into  the  identical  form  as  the  no  flow  case;  i.e.,  let 


1 


y2  = 1/(1  -M02) 

(K\  x',  y',  z')  = (yK,  yx.  y.  z) 

(0j.  <t>v  0)  = exp(-  i K'  M0  x'a  ) (0j,  0t.  0) 

cos  6 ' =(cos0-Mo)/(l  -M0cos  6) 
cos  0'Q  = (cos  dQ  + Mc)/< ! + Mq  cos  dQ  ) 

Substitution  of  the  transformed  quantities  into  equations  (C-3).  (C-4),  (C-7),  and  (C-8) 
results  in  the  following. 


0j  = exp  [i  K’  (- Xq  cosO'  + y-Q  sin0')| 

0t  = 0j  + V 


90* 

—=0  on  y’  = 0.  x'  < 0 
3y 

<V')2  0 +(K')2  0 =0 


(C-10) 


(C-ll) 


(C-12) 


where  (0j,  i//t,  0)  are  the  potential  functions  in  the  case  with  flow  for  the  incident,  total 
and  diffracted  fields,  respectively. 

These  are  identical  in  form  to  equations  (C-2),  (C-3),  (C-4),  and  (C-5)  for  the  no  flow 
problem.  Hence,  the  solution  for  (0t/0()  *8  given  by  the  static  potential  solution  (dq/dq) 
defined  with  transformed  coordinates: 


0t/0j  = 0t^i  I = G (r'  I ) 
transform 
coordinates 


(C-13) 


The  interpretation  of  this  solution  in  terms  of  pressure  or  velocity  potential  requires 
consideration  of  continuity  requirements  at  the  edge.  For  the  trailing  edge  (i.e.,  M0  is 
positive),  the  solution  (d»t/0i)  represents  the  pressure  field.  For  the  leading  edge  (i.e., 
M0  is  negative),  the  solution  represents  the  velocity  potential  field.  In  order  to  get 
pressure,  we  must  differentiate  the  solution  with  respect  to  the  space  variable  in  the 
direction  of  flow;  i.e., 


P i K'  0 - M0"^~7  0 
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Then  for  the  leading  edge 


APPENDIX  D 

THEORETICAL  DATA  CURVES 
FOR  CYLINDER  DIFFRACTION 


D.  G.  Dunn 


This  appendix  contains  the  theoretical  trend  figures  for  noise  diffracting  about  a 
fuselage  (cylinder)  structure.  A cross  reference  table  shown  below  is  provided  to  assist 
use  of  the  trend  figures  for  the  analysis  (sec.  2.2:2.2)  and  the  interpolation  procedure 
(sec.  4.3.1)  of  the  main  engineering  report,  volume  I. 


APPENDIX  E 

SOUND  REFRACTION  BY  AIRFOIL  TRAILING  WAKE 

H.  Y.  Lu 


This  appendix  covers  the  modeling,  the  formulation,  and  the  solution  of  the  problem  of 
sound  refraction  by  an  airfoil  trailing  wake. 


1.0  ANALYSIS 
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Sound  propagation  through  the  wing  trailing  wake  from  an  over-the-wing  engine 
presents  a problem  of  sound  reflection  and  refraction.  Theoretical  work  related  to  this 
problem  has  been  completed  at  The  Boeing  Company1  and  is  applicable  with 
modifications  to  the  analysis  of  sound  propagation  through  the  wing  trailing  wake. 

A few  simplifications  have  to  be  made  in  modeling  the  wake  and  the  noise  source.  The 
following  assumptions  are  made: 

1.  The  noise  sources  are  fixed  relative  to  the  airplane  and  can  be  represented  by  point 
acoustic  sources. 

2.  The  wing  trailing  wake  is  constant  in  velocity  and  temperature.  The  velocity  is 
different  from  that  of  the  freestream,  while  the  temperature  is  the  same. 

3.  The  wake  extends  to  infinity  with  finite  thickness  and  the  presence  of  the  w»ng  as 
a solid  body  is  ignored  in  this  part  of  the  analysis. 

4.  Results  from  this  analysis  can  be  superimposed  on  that  of  the  wing  diffraction 

5.  A perfect  and  inviscid  gas  is  assumed. 

These  assumptions  are  used  to  analyze  the  propagation  of  sound  through  the  wing 
trailing  wake. 

We  introduce  a time  dependent  point  source  Q 8(x)  8(y)  8(z)  in  the  flow.  Q has  a physical 
dimension  of  mass  per  unit  time.  The  continuity  equation  for  acoustic  distrubances  can 
be  written  as 


dS  ndS  du  dv  dw 
d t u dx  dx  dy  dz 


q 6(x)  6(y)  6(z) 


(E-l) 


‘H.  Y.  Lu,  "Acoustic  Far  Field  of  a Point  Source  in  Cylindrical  and  Parallel  Flow  Fluid  Layers." 
A1AA  Paper  No.  75-500,  March  1975. 
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where  q = p0,  and  that  S is  the  condensation  defined  by  p = p0(l  + S)  with  p0  as  the 

undisturbed  density  (fig.  E-l).  U is  the  velocity  of  the  uniform  flow,  u,  v,  and  w are  the 
components  of  acoustic  perturbation  on  velocity  along  x,  y,  and  z directions, 

Q 

respectively,  c is  the  speed  of  sound;  i.e.,  c = yPJpQ. 

The  momentum  equations  along  x,  y,  z directions  are 


dn  , , du 

at  dx 


= - c* 


2 dS 
dx 


iv  u^l 

dt  dx 


t as 
- c"  ay 


(E-2) 

(E-3) 


dw  . , dw  ,2  dS 

dt  U dx  C dz 


(E-4) 


The  dimensionless  acoustic  pressure  p is 


P-  P 
r ro 

P 


= y s 


(E-5) 


with  T = cp/cv 

From  equations  (E-l)  through  (E-4),  u,  v,  and  w are  eliminated  and  S is  substituted  by  p 
from  equation  (E-5).  The  result  is 

(ir+  u^)2 p - °2, v2  p = 7 (£+  u_lr)q‘ 6(x)  6(y)  6{z)  <e-6> 

where  V-  = -~r  + -^  + 

dx“  dy“  dz“ 


i 
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Apply  Fourier  transform  to  equation  (E-6), 

F"  + f2F  = ^2^(f'Mkx)  Q“6U) 


(E-7) 


where 


F(z) 


oo  oo  _ oo 

1 iff  -i(k„x+k  y-cot) 

1 / J J pe  y dxdydt 

8jt  _oo  Joo  — oo 


■ fff-e  i<kxx+kyy~<‘,t)  dkx  dky  dco 


P 


and 


2 2Mu>k  -)  * 

t2  = ~2  7 ‘X  -k2 

with  Mach  number  M = U/c.  Omega  (<u)  is  the  angular  frequency;  kx  and  ky  are  the 
wave  number  in  x and  y directions,  respectively. 

The  solution  of  equation  (E-7)  is 

if;Z  -if.-z 

Fj  = 3j  e i + bj  e J (E-8) 

where  j = 0,  1,  2,  3 for  regions  0,  1,  2,  and  3,  respectively. 

Since  regions  0,  1,  and  3 have  identical  flow  conditions, 

t i o , ,2  2M„  cj  kv  7 

?o  =H  = $3  --0-M5)k2  -k2  • 


cj~  2M2cokx 


1 

^ \ 1,* 


= O' - ( 1 - M2  >k£  - ky 


Eight  boundary  (Bi)  conditions  are  needed  to  determine  the  eight  unknown  constants  aj 
and  bj. 

Bl.  In  region  0,  we  have  only  the  outgoing  waves  and  therefore  a0  = 0 for  z 
B2.  Outgoing  waves  for  region  3,  l>3  = 0 as  z-*®. 

B3.  Integrating  equation  (E-7)  across  z = 0, 

F'](0)-F;(0)  = j1-X.(W-kxUo)Qa; 

4tt  C" 

or  more  explicitly, 

■fo  ^al  " bj ) + if0  bQ  = ^-(w  - kxUQ) 

4jr*-  c “ 


B4.  Continuous  pressure  at  z = 0, 

bQ=  a,  +b, 


Continuous  pressure  at  interfaces  z = zj,  and  z = zz- 
B5  a,/°''tb,c'ifoZU3  e‘f2^l  +b 
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B6  *?2  /2  . u '^2  Z1  _ z2 

DD-  a->  e + b ■>  e -a^e  u *■ 


Continuous  displacement  17  of  the  interface  at  z = Zj  and  z = z^,  we  have 


From  equation  (E-3) 


'd  a Y _ c2 

rt  + vTx)  '--^r 


After  Fourier  transform 


y (-  oj  + kxU)-  az 

where  t)  is  the  Fourier  transformed  displacement  which  has  to  be  equal  on  both  sides  of 
the  interface;  i.e., 

B7. 

, [a,  f„  .*•  *•-»,*  a,  - b,  f, 

y (co  ~ kY  UJ2  L 1 0 J (w  - k„  U-,)-  L “ " J 


C2  T V ‘f-*  z2  . -ifr  Z?1  c2  , ifo z? 

,",V  a2  h e -b2f2e  “ = ~ TTT7  a3foe 

- k..  LM~  L J 7(u-kxU0)- 


7 (co  - kx  U2) 


These  eight  boundary  conditions  provide  the  solution  of  the  coefficients  at  and  83  given 
at  the  end  of  this  appendix.  In  region  3 we  have 

p *o  2 

F3  = a3e 

The  dimensionless  pressure  for  each  frequency  a>  is 


ff 


a3  e dkx  dk\ 


where  we  have  made  the  following  change  of  variables  (fig.  E-l). 
x = R sin  if  cos  9 

y = R sin  y sin  9 

z = R cos 


-[(f)2- ^ 

+ ky  sin  <f  sin  9 


- ( 1 - Mq4-)  kx^  - ky“  cos  ifi  + kx  sin  <p  cos  9 


The  integrations  in  equation  (E-9)  can  be  performed  asymptotically  for  R-b®. 
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2.0  FAR-FIELD  SOLUTION 
AND  THE  TRANSMISSION  COEFFICIENT 


When  the  observer  is  at  a great  number  of  wave  lengths  away  from  the  source,  an 
asymptotic  solution  for  the  acoustic  far/field  can  be  found  by  the  stationary  phase 
method.  Equation  (E-9)  is  evaluated  at  the  stationary  phase  for  the  asymptotic 
expression  as  R-»°o 

1 (1  - M2  ) sin2ip  cos2d 

_ 2xri  co  , I_Mo  !_m2  +Mo  sin2*>cos20  iR  h3c 

P3co  ~ R r ^ ■■  — 5 a3c  c ~ (E-10) 

1 - sin“V  cos ‘■d 

where  83s  and  h3s  are  a3  and  I13  evaluated  at  the  stationary  phase.  The  stationary 
phase  wave  numbers  are  given  in  the  following  section. 


We  define  the  sound  transmission  coefficient  Tco  for  the  frequency  co  as 

T ~ P3co 
<0  D 

Pico00 


(E-ll) 


where  Pl^oo  is  the  reference  dimensionless  acoustic  pressure  for  the  same  source  in  an 
infinite  air  at  rest.  We  have 


2ni  u>  i R co  /c 

P | ojoo  ~~  R c ^ 00 


(E-12) 


where 


^loo  87^*-  C cos 


For  a subsonic  flow  case  we  have 

( 1 _ \ cin2 


, (1  - M2  ) sin2^  cos2d 

1 - M2  1 - M2  + M2  sin2v?cos2d  a3S 

i T”  ~ Z 

1 - sin*v  cos^o  d loo 


(E-13) 


For  convenience,  let  o>  = 2wf,  where  f is  the  frequency.  The  transmission  coefficient 
takes  the  form 


If  = 20i°8|0  1 T2it  f 


which  gives  a unit  of  decibels. 


urn  ■ 


t 


3.0  NOTES  ON  EXPRESSIONS  FOR  a1(  a3 

AND  STATIONARY  PHASE  WAVE  NUMBERS 


1 7 (<^-kx  U0) 


al  ~ ? ~7 


a3 


8tt 
' 0 


c-  r 


to 


C 2 ^ 


4 n2  7(oj-kx  UQ)(cj*kx  U2)2  A sin(f2D^  + iB  cos^2D) 


where 


D = z2  - z i 

cV, 


c4f  2 

V * O 


7"-(w-kxU2r  7"  (co  - kx  UQr 


2C4f2f0 


72(w-kx  U2)2(co-kx  UQ)2 


k = — 

KXS  r 


sin  <p  cos  9 


L(l-M2  + M2  sin2*  cos2*)**  1-M2J 


k = -^ 
ys  C 


y -y  y 

I ( 1 - ) sin“<£  cos“0 


1_1  - M“  I - M2  + M2  sin2t0cos20  J 


'A 


sin  ■p  sin  6 


-)  i Vi 
(1  - sin “10  cos -0) 
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APPENDIX  F 

SPATIAL  AVERAGING  FACTOR  I 

C.  H.  Berman 


I 


i 

I 

I 

I 

1 


1 

i 


< 


j 


11 


The  averaging  factor  I is  defined  as 


u'“  dqs  dx  dy 

UA^  <x  + + y~  + h" 


id  ds 

In  2 dft 

Lua 


typical 


with  the  source  located  at  x = -b,  y = h,  z = 0 relative  to  the  trailing  edge  as  shown  in 
figure  37  of  volume  I. 


We  will  arbitrarily  write 


typical 


1 


c 


and  consider  d<rs/dfl  which  is  actually  a function  of  8=  8(x)  at  fixed  frequency,  to  be 
nearly  constant  as  a function  of  x.  This  is  almost  true  for  the  higher  Strouhal  number 
(S)  values  shown  in  figures  41  and  42  of  section  2. 2.4. 2 (vol.  I),  and  our  estimates  will 
involve  cases  which  are  either  in  or  close  to  this  regime.  Direct  integration  then  leads 
to 


For  example,  in  the  limit  b = h = c 


I = fin 


v^T+  1 


0.881 


for  h = 0 


c 


which  has  the  value  1=1  when  b = c and  a logarithmic  singularity  when  b/c-+0. 


i 
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Thus,  scattering  evaluations  using  I * 1.0  will  be  fairly  good  unless  the  source  starts  to 
get  close  to  the  trailing  edge  where  a logarithmic  singularity  occurs. 


'! 

i 

\ 
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APPENDIX  G 

SPECTRUM  SHAPE  FACTOR 
FOR  JET/EDGE  INTERACTION  NOISE 


r \ 


L.  Filler 


The  normalized  sound  pressure  level  given  by  equation  (82)  of  volume  I is  defined  to 
have  a maximum  at  f/fQ  = 1;  i.e.,  f0  is  the  frequency  at  which  the  spectrum  peaks. 
Therefore, 


'LS(f/f0)] 


d(f/fQ) 


0,  at  f/f0  = 1 


Differentiating  equation  (82)  (vol.I)  gives, 

d[s(f/f0)]  [l  + F(f/f0)7]‘' 


d(f/f0) 


= 10 


( i + F)va/f0p 


(1  + 


d(f/f0) 


"V 


l[l  +F*(l7f0i7_ 


= 0 


which  is  satisfied  by, 


(l~J  d [l  +F(f/f0)7]1'  - [l  +F(f/f0>7]\lj(f7 


l7f  0)M  }=  0 


Carrying  out  the  indicated  differentiations  yields,  after  a little  algebra, 

(7  v-n)  F (f/f0)7  - n = 0 


For  the  maximum  to  occur  at  f/f„  — 1,  the  required  condition  is, 


F = 


7 V-n 


(G-l) 


It  is  immediately  seen  that  F is  the  ratio  of  the  low  frequency  (f/f0  <<  1)  to  high 
frequency  (f/f„  >>  li  falloff  from  the  peak.  From  equation  (82)  (vol.  I),  the  low  frequency 
limit  is. 


Sit  l0l  * 10  lo>i  I * I i'  (t  t0)^  — < < 


(G-2) 


and  the  high  frequency  limit  is 

»»  «„»  * " i ‘ ) 

hence,  fi-Hyi'  i is  the  staled  rati.. 


^ > '• 


(G-3) 


I 
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APPENDIX  H 

RELATIVE  VELOCITY  SCALING 
IN  JET/EDGE  INTERACTION  NOISE 


L.  Filler 


I! 


i 

1 

i 


For  simplicity,  consider  forward  velocity  effects  for  the  nozzle  on  the  wing  case.  From 
the  estimation  equation,  the  relevant  terms  are 

1 0 ]og(v j v£  )-IO 
and  two  cases  may  be  distinguished: 


1.  If,  L « LCo  it  follows  that  L =s  LCo  Vj/Vr  is  always  satisfied.  In  this  case  the  abov  • 
terms  can  be  written  as, 


10 


10  log  V 


5 

R 


which  confirms  the  result  of  Bhat  and  Gallo-Rosso  (ref.  49,  vol.  I)  that  for  L s L(. 
the  interaction  noise  is  reduced  with  forward  velocity  and  depends  on  relative 
velocity  to  the  fifth  power. 


2.  If,  L > LCo  Vj/VR  there  results 


lOlog^Vj  VR4) 


10  log  V?  VR 


which  states  that  the  interaction  noise  does  not  decrease  as  rapidly  but  depends  on 
the  product  of  jet  velocity  squared  and  relative  velocity  cubed  when  the  stretched 
potential  core  length  remains  smaller  than  the  shield  length. 
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